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1.0  HITRAN  DATABASE 


Work  done  under  this  contract  produced  several  major 
changes  to  the  AFGL  line  parameter  compilations.  First,  the  two  line 
parameter  compilations  (the  main  and  trace  gas  atlases)  were 
merged  into  a  single  entity  and  called  high-resolution  transmission 
molecular  absorption  database  (HITRAN).  Next,  the  format  of  the 
database  was  changed  and  other  parameters  corresponding  to  each 
quantum  transition  were  added.  Each  record  of  a  transition  now 
consist  of  100  characters  of  information  rather  than  the  previous  80 
characters.  Finally,  the  database  is  now  distributed  with  several 
other  files  that  assist  in  reading  and  using  the  database  and  which 
give  useful  information  relating  to  the  database. 

The  first  version  of  the  combined  databases  was  issued  in 
late  1986  and  was  called  the  HITRAN  database:  1986  edition.  It  was 
formed  from  merging  the  1982  editions  of  the  AFGL  main  gas 
compilation1  and  the  trace  gas  compilations. 2  In  merging  these  two 
compilations,  much  new  data  was  added  with  the  corresponding 
previous  data  being  removed.  All  of  these  changes  were  described  in 
an  accompanying  publication. 3  Originally  the  database  contained  the 
following  molecular  parameters  for  each  quantum  transition:  (1) 
resonant  wavenumber;  (2)  line  intensity;  (3)  air-broadened 
halfwidth;  and  (4)  lower  state  energy  (as  well  as  the  unique  quantum 
state  identifications).  New  parameters  have  been  added  to  this  list  to 
allow  the  database  to  be  applicable  to  remote  sensing  in  the 
atmosphere  and  capabilities  to  deal  with  non-local  thermodynamic 
equilibrium  effects  in  the  upper  atmosphere.  The  new  parameters 
are:  the  transition  probability,  the  self-broadened  halfwidth,  and  the 

exponent  of  temperature  dependence  of  the  air-broadened 
halfwidth.  Data  fields  have  also  been  reserved  for  the  pressure  shift 
of  the  transition,  accuracy  criteria  for  the  three  principal  parameters, 
and  references  to  the  sources  of  the  latter  parameters.  The  card 
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image  concept  (80  characters/transition)  was  abandoned  and  100 
characters/transition  was  adopted.  In  order  not  to  exceed  100 
characters,  the  isotope  labels  and  vibrational  state  labels  were 
indexed.  Recalling  of  the  information  is  accomplished  via  the 
program  SELECT  which  is  distributed  with  the  database.  This 
program  is  discussed  later. 

The  creation  of  a  new  database  gave  us  an  opportunity  to 
clean  up  many  of  the  minor  errors  that  were  present  on  the  1982 
versions.  Many  of  these  errors  were  not  serious,  e.g.  many  species 
had  the  quantum  designations  present  in  several  formats,  different 
symmetry  labels  were  used  for  the  same  designations,  etc.  Some 
errors  were  more  serious,  parts  of  or  the  entire  quantum  field  was 
missing,  thus  the  data  can  not  be  associated  with  a  specific  transition. 
The  reformatting  allowed  us  to  define  six  groups  of  rotational 
quantum  number  formats  for  the  molecules  and  ensures  that  the 
1986  version  will  have  a  single  format  for  each  species.  The 
rotational  formats  are  given  in  Table  1. 

One  of  the  new  parameters  being  added  is  the  transition 
probability,  Rjf,  in  Debyes  squared.  This  parameter  is  independent  of 
both  temperature  and  isotopic  abundance.  The  transition  probability 
is  related  to  the  intensity  of  a  transition,  Sj f ,  from  state  i  to  state  f  by 


Sj  KT)  =  V: 


3hc 


if 


1  -  exp( 


-C-)V 


5^) 

/T  . 


gjl  a 

Q(T) 


exp(*c2E^r)  Rjf*  10 


-36 


(1) 


Here  vjf  is  the  resonant  frequency  of  the  line  in  wavenumbers,  Ej  is 
the  lower  state  energy  of  the  transition,  gj  is  the  nuclear  spin 
degeneracy  of  the  lower  level,  Q(T)  is  the  total  internal  partition  sum, 
Ia  is  the  natural  isotopic  abundance,  and  C2  is  the  second  radiation 
constant  (hc/k).  The  transition  probability  on  the  database  was 
obtained  by  removing  all  terms  on  the  right  hand  side  of  eq.  (1) 


Group  1:  Asymmetric  Rotors, 


H;0,  O3, 

S0?  , 

no2  , 

HN03, 

h2 

CO, 

H0C1 , 

h2o2. 

h2s 

J'.  Ka', 

V  • 

F'  . 

Svm'  ; 

J" 

.  Ka 

ir  prt 

Svm 

12.  12  . 

12  , 

12, 

Al  : 

12 

,  12 

,  12 

,12, 

Al 

Group  2:  Diatomic  and  Linear  Molecules  with  Integer  J 

CO,.  N90,  CO,  HF,  HC1,  HBr ,  HI,  OCS ,  N2 ,  HCN ,  C2H2 

Br.  F" ,  _ Br,  J",  Sym" 

5X , Al ,  1 2  , IX ;  4X , Al ,  13,  Al 

Group  i:  Spherical  Rotors 

Methane  (CH,t  only,  not  CH3D) 

J',  R'  ,  C',  N:'  ,  Sym';  J",  R" ,  C" .  N" ,  Svm" 

12,  12.  A2,  12,  Al  ;  12,  12,  A2,  12.  Al 

Group  *:  Symmetric  Rotors 

CH.D,  Nil,,  CH3CI,  C2H6*,  PH 3 

J'.  K',  C',  ,  Svin'  ;  J",  K"  ,  C"  ,  ,  Sym" 

12,  12.  A 2 ,  2 X ,  Al  ;  12,  12,  A2,  2X,  Al 

Group  0:  Triplet  (.'round  Electronic  States 
°2 

Br,  F"  ,  Br,  N’ " ,  Br,  J”,  Sym" 

3  X .  A 1 ,  F  4 . 1 ,  IX;  Al,  12,  Al,  12,  2X,  Al 

Group  6:  Doublet  Ground  Electronic  States  (Half  Integer  J) 

NO,  OH.  CIO 

.  Br,  F”  ,  ;  ,  Br ,  J”  .  Sym" 

SX ,  Al.  12,  IX:  3X ,  Al,  F4.1.  Al 

Prime  and  double  primes  refer  to  upper  and  lower  states 
respectively;  Br  is  the  P- ,  Q- ,  or  R-  branch  symbol:  J  is  the 
rota; ional  quantum  number;  Svm  is  e  or  f  for  1-type  doubling,  +  or 
-  tor  '.rmme ; rv ,  etc.  (for  further  explanation  see  references). 

stnt.lv  using  group  2.  quantum  notation  on  the  tape. 

- u L 1"  : .  format  of  Rotational  Quantum  Numbers 
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(except  Q(T))  out  of  Sjf.  (The  partition  sums  were  not  ail  known  at 
the  time  of  producing  this  version.)  The  data  for  the  line  strength, 
resonant  frequency,  and  lower  state  energy  come  off  the  database. 
The  reference  temperature  of  the  database  is  296  K.  The  isotopic 
abundances  used,  now  supplied  as  part  of  the  third  file  distributed 

with  the  database,  are  given  in  Table  2.  The  nuclear  spin  degeneracy 

factor,  gj,  depends  on  the  particular  isotope  symmetry.  These  were 
calculated  for  each  isotopic  variant  on  the  database. 

When  there  are  equivalent  nuclei  in  a  molecule,  i.e. 

identical  isotopes  of  the  same  element  which  have  exactly  the  same 
molecular  environment,  the  nuclear  spin  statistics  will  affect  the 
populations  of  the  molecular  states  and  hence  the  transition 

intensities. 

Consider  first  two  equivalent  nuclei.  The  molecule  has  a 
two-fold  axis  of  symmetry,  and  coordinates  of  the  two  equivalent 
nuclei  may  be  interchanged  by  a  180°  rotation  about  this  symmetry 
axis,  or  by  combinations  of  inversion  and  180°  rotations  about 

various  axes.  If  the  nuclei  obey  Bose-Einstein  (integer  spin) 
statistics,  an  interchange  of  both  the  position  and  spin  coordinates 

must  leave  the  wave  function  unchanged;  if  they  obey  Fermi-Dirac 
(half-integer  spin)  statistics,  the  wave  function  must  change  sign.  It 
can  be  shown^  that  the  total  number  of  symmetric  spin  functions 
possible  (hence  the  weight)  is 

nsym  =  8s  =  (21  +  1)0  +  1)  (2a) 

and  of  antisymmetric  spin  functions  possible  is 

nantisym  =  8a  =  (21  +  1)  (I).  (2b) 

A  molecular  rotation  of  180°  about  the  symmetry  axis  and  an 
exchange  of  spin  coordinates  of  the  two  equ  valent  nuclei  amounts  to 
an  exchange  of  all  coordinates  of  these  two  nuclei.  Hence  if  Bose- 

Einstein  statistics  apply  to  the  two  nuclei,  the  symmetric  spin 


r  ~  1  n  L 

i  v  e  M  a  ‘  .i  i 

1  a  1  A 3  1 1  e  f  * " 1  ]  ■  i i  '  .'i  l  F  p  e 

c  ic-s  on  Database 

1  ••rule 

A  f  G  L 

lelative  "cl e  rule 

AFGL  relative 

isotope 

na  t  u  t  a  1 

isotope  natural 

code 

a  b  u  n  d  a  n  c  e 

code  abundance 

H20 

161 

0.9973 

H  NO  3 

146 

0 .9891 

181 

0 .0020 

1  71 

0 .0004 

OH 

61 

0.9975 

162 

0.0003 

81 

0.0020 

62 

0.00015 

CO  2 

626 

0 . 9842 

636 

0.0110 

HF 

19 

0.99985 

628 

0.0039 

627 

0 .0008 

HC1 

15 

0.7576 

638 

0 .000044 

17 

0.2423 

637 

0 . 0000G9 

828 

0.0000040 

HR  r 

19 

0 . 5068 

/  2  8 

0 .000002 

11 

0.4930 

03 

666 

0.9928 

HI 

17 

0.99985 

668 

C  .  00  4  0 

686 

0 . 0020 

CIO 

56 

0.7559 

76 

0.2417 

N20 

4  46 

0.9904 

4  36 

0.0036 

OCS 

622 

0.937 

5  46 

0.0036 

624 

0.0416 

4  48 

0 .0020 

632 

0.0105 

447 

0 .0004 

822 

0.0019 

CO 

26 

0 .9865 

H  2  CO 

126 

0.9862 

36 

0.011 

136 

0.0111 

28 

0  .  0020 

128 

0.0020 

27 

0 .0004 

HOC  1 

165 

0.7558 

CH  4 

211 

0.9883 

167 

0.2417 

311 

0.0111 

212 

0 .00059 

N2 

44 

0.9928 

02 

66 

0.9952 

HCN 

124 

0.9852 

6  8 

0.0040 

134 

0.0111 

67 

0 .0008 

125 

0.0036 

MO 

46 

0.9940 

CH3C1 

215 

0.7490 

56 

0.0036 

217 

0.2395 

48 

0 .0020 

H  202 

1661 

0.9949 

S02 

626 

0.9454 

646 

0.0420 

C  2  H  2 

1221 

0.9776 

1231 

0.0219 

::o2 

646 

0.9916 

C2H6 

1221 

0.9776 

:;h  3 

4  111 

0 . 9960 

5111 

0 . 0036 

PH  3 

1111 

0.99955 

Table  2 


functions  must  be  used  with  rotational  functions  which  are 

symmetric  with  respect  to  a  rotation  around  the  symmetry  axis,  the 
antisymmetric  spin  functions  with  antisymmetric  rotational 
functions.  The  opposite  is  the  case  for  Fermi-Dirac  statistics.  The 
ratio  of  the  number  of  spin  states  or  the  statistical  weights  of  the 
levels  which  are  symmetric  to  those  which  are  antisymmetric  with 

respect  to  a  180°  rotation  about  the  symmetry  axis  is, 

I  +  1 

For  Bose-Einstein  statistics:  1 

1 

For  Fermi-Dirac  statistics:  1+1.  (3) 

The  line  strength  for  molecules  with  two  equivalent 
atoms  contains  the  nuclear  spin  statistical  weight  factor  from  eq.  (2). 

If  there  are  more  than  one  pair  of  equivalent  nuclei  in  a 

molecule  the  symmetry  properties  of  each  pair  may  be  taken  into 

account.  An  example  is  CH2CI2.  If  we  let  the  spins  of  the  two 
different  types  of  nuclei  be  1 1 ,  and  I2  we  have 

nsym  =  (21 1  +  I)(2l2  +  1  )(21 1 12  +  Ij  +  I2  +0  (4a) 

and 

nanti  =  (21 1  +  1)(2I2  +I)(21il2  +  I]  +  12)-  (4b) 

In  general,  for  a  molecule  with  n  pairs  of  identical  nuclei,  the  number 
of  symmetric  spin  functions  which  can  be  formed  is 

fl(2Ij+l)+1 

J  (5a) 

and  the  number  of  antisymmetric  functions  is 


f 


nanti  —  -) 


(5b) 


fl|2>,+  1) 

npi.+'H 

—> 

II 

i=i 

If  the  molecule  has  a  two-fold  axis  of  symmetry  and  the 
n  pairs  of  nuclei  are  interchanged  by  a  180°  rotation  around  this  axis, 
the  ratio  of  intensities  of  rotational  levels  of  odd  and  even  K  is  given 
by  eq.  (5). 

The  cases  described  above,  i.e.  180°  rotations,  apply  to 

linear  and  asymmetric-top  molecules.  Next  we  consider  symmetric 

top  molecules  which  may  have  three  or  four  equivalent  nuclei  in  a 
molecule.  For  molecules  with  three  equivalent  nuclei,  the  levels  with 
K  a  multiple  of  three  always  have  greater  weight.  The  statistical 
weights  for  these  levels  are  given  in  Table  3.  Similarly  for 

symmetric  top  molecules  with  four  identical  nuclei  equidistant 

around  the  axis,  the  statistical  weights  are  in  Table  4. 

The  nuclear  statistical  weights  for  species  present  on  the 
atlas  can  now  be  given. 

H2O  This  is  an  asymmetric  rotor  with  two  identical  nuclei  for  the 
species  161,  171,  and  181.  From  eq.  (2)  we  have 

nsym  =  8s  =  (2I+1)(I+1);  nanjj  =  ga  =  (2I+1)(I). 

For  hydrogen  I  =  1/2,  which  gives  gs  =  1  and  ga  =  3  where 
the  even-(odd)ness  is  determined  by  the  even  or  odd-ness 
of  Ka  +  Kc  +  V3.  For  HDO  there  are  no  equivalent  nuclei,  thus 
8s  =  8a  =  1- 

CO2  There  are  eight  isotopic  species  in  the  atlas  for  CO2.  Of 
these,  only  three  have  a  pair  of  identical  nuclei,  626,  636, 
and  828.  All  other  species  do  not  have  interchangeable 
nuclei  and  ga  -  8s  =  *  For  the  species  which  have  nuclear 
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Table  3.  Statistical  Weights  for  Symmetric  Tops  with  Three 

Identical  Nuclei 


State  Statistical  Weight 

K  a  multiple  of  3,  but  not  0  1/3  (21  +  1  )(41  —  +41  +  3) 

K  not  a  multiple  of  3  1/3  (21  +  1  )(41—  +  41) 
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Table  4.  Statistical  Weights  for  Symmetric  Tops  with  Four  Identical  Nuclei 


K  =  0 


( 1  + 1  )(2I  + 1  )(2I2+I  + 1 )  for  K  a  multiple  of  four 


Bose-Einstein 
statistics  or 


I(2I+l)(2I2+3I+2) 


K  even,  for  K  not  a  mul¬ 
tiple  of  four 


Fermi-Dirac 

statistics 


for  K  a  multiple  of  four  Fermi-Dirac 

statistics 


K  even,  K  not.  a  multiple 
of  four 


Bose-Einstein 

statistics 


1(1  +  1  )(2I  +  1  )2  K  odd 


K  =  0 


Bose-Einstein  Statistics 
1/2  (1+1  )(2I  +  1  )(2I2  +  3I  +  2)  J  even 

1/2  I  (I+1)(2I-1  )2I+1)  J  odd 


Fermi-Dirac 
1/2  1  (21+1  )(2I2  +  I  + 1 ) 
1/2  I  (I+l)(2I+l)(2I+3) 


J  even 
I  odd 
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weights,  we  note  both  16q  and  ^0  have  I  =  0  so  all  three 
species  have  the  same  nuclear  statistical  weights,  namely 

gs  =  (2(0)+l)(0+l)  =  1  and  ga  =  (2(0)+l)(0)  =  0. 

All  CO2  lines  on  the  database  have  gj  =  1. 

O3  The  668  species  has  no  interchangeable  nuclei  thus  g  =  1, 

the  other  two  species  have  as  the  equivalent  nuclei 

related  by  a  180°  rotation  hence  from  eq.  (2)  gs  =  1  and  ga  = 
0.  For  ozone  all  lines  present  on  the  atlas  have  gj  =  1. 

N2O  This  molecule  has  the  structure,  NNO,  and  there  are  no 

nuclei  that  are  equivalent.  There  are  many  molecules  on  the 
atlas  that  are  similar  to  this  and  are  listed  below.  Note  for 
all  isotopes  of  these  species  there  are  no  nuclear  statistical 
weight  factors  to  consider  for  coupling  with  rotational  wave 
functions:  gj  =  1  for 

N20,  CO,  NO,  OH,  HF,  HC1,  HI,  HBr,  OCS,  HOC1,  and  HCN. 

CH4  Methane  is  a  symmetric  top  molecule  with  four  equivalent 
nuclei  for  the  atlas  species  211  and  311.  The  nuclear 

statistical  weights  can  be  obtained  from  Table  4  or  can  be 
related  to  the  labels  used  to  describe  the  state.  A,  E  or  F. 
This  gives  the  following  association  for  gj:  for  A  species  gj  = 
5,  for  E  species  gj  =  2  and  for  F  species  gj  =  3. 

The  deuterated  species,  CH3D,  has  three  equivalent  nuclei 
which  gives  from  Table  3  with  I  =  1/2 

gj  =  4  for  K  a  multiple  of  3 

gj  =  2  K  not  a  multiple  of  three  (includes  0) 
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O2  Here  the  two  oxygen  atoms  are  equivalent  for  the  66  species 
which  gives  gs  =  1  and  ga  =  0,  and  not  equivalent  for  the  67 
and  68  species  so  gj  =  1. 

SOj  The  two  species  present  on  the  atlas  626  and  646  have 

equivalent  1^0  nuclei  hence  gs  =  1  and  ga  =  0. 

NO2  The  646  species  is  the  only  one  present  on  the  atlas  and  has 
gs  =  1  and  ga  =  0. 

NH3  Ammonia  is  a  symmetric  top  with  three  equivalent  H  atoms 
for  both  species  on  the  atlas.  From  Table  3  we  find 

K  a  multiple  of  three  but  not  0  g  =  4 

K  not  a  multiple  of  three  (including  0)  g  =  2 

HNO3  The  structure  of  nitric  acid  is  such  that  there  are  no 

equivalent  l^O  nuclei.  This  implies  that  gj  =  1. 

H2CO  All  three  species  on  the  atlas  contain  a  pair  of  equivalent 

hydrogen  nuclei.  We  have  for  levels  with  Ka  +  Kc  +  V3  even 

gs  =  1  and  with  Ka  +  Kc  +  V3  odd  ga  =  3. 

N2  The  nitrogen  atoms  are  symmetric  by  180°  rotation,  1(1  ^N)  = 
1  so 

gs  =  (2(1)  +  1)(1  +  I)  =  6 
and 

ga  =  (2(1)  +  1)(1)  =  3. 

CH3CI  Both  species  on  the  atlas  have  three  identical  hydrogen 

nuclei,  thus  we  have 
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g  =  4  when  K  is  a  multiple  of  three  but  not  0  and 

g  =  2  when  K  is  not  a  multiple  of  three,  including  0. 

H2O2  In  this  cis-conformation  a  rotation  through  180°  exchanges 
both  the  oxygens  and  hydrogens.  In  the  trans-conformation 
inversion  through  the  center  of  gravity  of  the  molecule  does 
the  same.  The  nuclear  statistics  are  given  by  eq.  (4)  where 
II  =1/2  and  I2  =  0.  The  even-oddness  is  determined  by  Ka 
+  Kc  +  V3  for  the  state 

8s  or  even  =  3  and  ga  or  odd  =  1. 

C2H2  Acetylene  when  rotated  180°,  exchanges  the  carbon  and 

hydrogen  atoms  so  again  we  use  equation  (4)  where  Ij  = 
1/2  and  I2  =  0  to  give  gs  =  3  and  ga  =  1. 

^2^6  The  structure  of  the  molecule,  H3C  -  CH3,  is  such  that  a 

rotation  along  the  internuclear  axis  by  120°  or  240°  leaves 
the  molecule  unchanged.  Thus  the  hydrogen  atoms  at  either 
end  are  equivalent.  This  gives  g  =  4  when  K  is  a  multiple  of 
3  and  g  =  2  when  K  is  zero  or  not  a  multiple  of  three. 

PH3  This  has  three  identical  nuclei  with  spin  1/2  thus  g  =  4  when 

K  is  a  multiple  of  three  and  g  =  2  when  K  is  zero  or  not  a 

multiple  of  three. 

The  other  new  parameter  added  to  the  database  is  the 
temperature  dependence  of  the  air-broadened  halfwidth.  This  is 
called  the  temperature  exponent  because  it  appears  as  an  exponent 
of  the  temperature  ratio.  This  becomes  clear  once  the  formulation  is 
considered.  The  halfwidth  can  be  considered  as  a  function  of  density 
or  pressure,  i.e. 

y  proportional  to  p  v  o  (density)  (6a) 


or 


1? 


y  proportional  to  p  v  o  (pressure). 


(6b) 


In  the  first  case  the  halfwidth  has  the  units  cm*1  /atm  or  MHz/torr, 
the  second  being  cnr1  or  MHz.  Formula  (6a)  leads  to  the  following 
temperature  dependence  (assuming  the  optical  cross-section,  a,  is 
temperature  independent) 


To  remove  the  negative  sign  of  the  exponent  this  is  often  rewritten 
as 


When  the  temperature  dependence  of  the  optical  cross-section  is 
added  (assuming  the  form  a(T)  =  Tm)  the  resulting  formula  is 


Working  in  the  other  formulation  (eq.  6b)  gives  (o  temperature 
independent) 


The  tilde  above  y  is  to  indicate  the  pressure  dependent  halfwidth  in 
cm-'  or  MHz  units.  The  confusion  arises  in  the  similarity  between 
eqs.  (9)  and  (7a)  in  the  exponents.  Only  on  closer  inspection  does  one 
find  the  temperatures  are  inverted.  When  the  temperature 
dependence  of  c  is  added  to  eq.  (9)  we  write  the  exponent  as  m  and 
note  that  m  =  n  +  1.  The  formulation  most  often  used  is  eq.  (7a)  with 
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the  halfwidths  given  in  (frequency  per  unit  pressure)  units.  Here  we 
use  equation  (8)  to  describe  the  temperature  dependence  so  we  need 
not  worry  about  the  negative  sign.  With  this  in  mind  we  now  review 
the  temperature  exponents  to  be  used  on  the  atlas. 

Table  5  lists  the  temperature  exponents  used  on  the 
database.  We  compare  our  values  with  those  found  on  the  GEISA 
atlas. 5  The  values  used  on  the  database  are  from  references  6-20. 
For  water,  GEISA  states  a  value  of  n  =  0.65  for  pure  rotation  from  the 
work  of  Bauer  et  al^  and  n  =  0.5  for  all  other  bands  from  the  work  of 
Davies  and  OIL?  In  reference  6  we  found  that  the  average  value  of  n 
=  0.65  was  for  self-broadening  of  H2O  from  ATC  and  QFT  calculations 
cited  in  the  references.  Bauer's  experimental  measurements  for  self¬ 
broadening  given  n  =  1.17  ±  0.06.  For  air-broadening  Bauer  reports 
"the  size  of  the  errors  is  such  that  the  temperature  dependence 
cannot  be  determined"  and  no  value  is  given.  The  reference  of  n  = 
0.5  from  Davies  and  Oli  is  merely  a  statement  of  the  classical  value  of 
n  (i.e.  the  assumption  of  a  temperature  independent  optical  cross- 
section).  Thus  both  reported  values  are  not  correct. 

In  the  Davies  and  Oli  paper,  four  transitions  were  studied 
for  temperature  dependence,  three  pure  rotation  and  one  V2  line. 
The  average  exponent  for  N2 -broadening  of  the  pure  rotation  lines 
was  n  =  0.64,  and  the  value  for  the  V2  line  was  n  =  0.45.  We  have 
discussed  the  results  with  Dr.  Richard  W.  Davies  and  feel  that  the 
value  for  the  642— >  753  V2  transition  (n  =  0.45)  represents  more 
how  the  temperature  exponent  can  vary  from  line  to  line  rather  than 
being  indicative  of  the  average  value.  Until  further  calculations  are 
made  it  is  suggested  the  value  n  =  0.64  be  used  for  all  water  lines. 

The  compilation  will  have  several  other  temperature 
exponents  that  do  not  agree  with  the  GEISA  values.  This  arises  from 
recent  measurements  not  yet  added  to  the  GEISA  database. 


Molecule 

HVHM 

min. 

max . 

n 

min.  max. 

H,Q 

0.0061 

0.1046 

--  0.64  -- 

CO, 

0.0559 

0.0899 

0.75  0.79 

03 

0.045 

0.077 

--  0.76  -- 

N,0 

0.0636 

0.0974 

0.64  0.82 

CO 

0.047 

0.088 

--  0.69  -- 

ch4 

0.0445 

0.0920 

0.63  1.00 

0, 

0.032 

0.062 

--  0.5  --- 

NO 

0.043 

0.063 

--  0.5  --- 

SO, 

0.110 

0.152 

--  0.5  --- 

0.062 

0.073 

--  0.5  --- 

nh3 

O 

O 

0.090 

--  0.5  --- 

KNO  3 

--  0.13 

--  0.5  --- 

OH 

--  0.083 

--  0.5  --- 

HF 

0.020 

0.126 

--  0.5  --- 

HC1 

0.0135 

0.0953 

0.20  0.88 

HBr 

0.015 

0.123 

--  0.5  --- 

HI 

0.008 

0.100 

--  0.5  --- 

CIO 

--  0.085 

--  0.5  --- 

ocs 

--  0.07 

--  0.5  --- 

h2co 

0.107 

0.108 

--  0.5  --- 

HOC! 

--  0.06 

0.5  --- 

n2 

--  0.06 

--  0.5  --- 

HCN 

0.0819 

0.1566 

--  0.5  --- 

CHjCl 

--  0.08 

--  0.5  --- 

H202 

--  0.10 

--  0.5  --- 

c,h2 

0.0400 

0.1158 

--  0.75  -- 

c2h. 

--  0.10 

--  0.5  --- 

ph3 

--  0.075 

--  0.5  --- 

R'ir.~.e  of  Air-Fr' -i  ioned  Ha If widths  and  Temperatur 
L  oo.'T.dencos 
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For  ozone  the  recent  work  of  Gamache^  gives  a  value  of 
n  =  0.76  for  air-broadening  temperature  dependence.  For  CO,  there 
are  several  references;  the  work  of  Bouanich  et  alJ6  giving  N  =  0.75 
and  n  =  0.77  for  the  temperature  ranges  200  -  300  K  and  93  -  300  K. 
The  work  of  Hartmann  et  al.^  giving  njq2  =  0-69  ±  0.02  and  that  or 
Connor  and  Radford^  giving  najr  =  0.74  ±  0.10.  We  have  decided  to 
use  the  value  reported  by  Hartmann  et  al.  for  the  temperature 
exponent,  n  =  0.69. 

For  NO2,  measurements  by  Devi  et  al.,9-20  for  ten 
transitions  give  an  average  value  of  n  =  0.97.  This  will  be  used  for 
all  NO2  lines. 

With  the  new  parameters  defined,  we  reformatted  the 
two  atlases  and  merged  them  into  one  database.  This  procedure 

employed  the  use  of  many  programs  to  check  for  the  consistency  ot 
the  data,  proper  formats,  etc.  From  the  many  checks,  we  are 
reasonably  confident  the  database  is  free  from  format  errors. 
Because  of  the  new  structure  of  the  database  and  the  adoitional 

parameters,  the  1986  version  is  distributed  with  a  user  friendly 
program  that  will  extract  information  from  the  database  and 
generate  files  of  several  specific  forms.  The  intent  was  to  ease  the 
change  to  the  new  form  for  users  conditioned  to  the  old  format.  The 
program,  called  SELECT,  uses  help  screens  in  the  interactive  mode  to 
prompt  the  user  for  specifics.  Any  incorrect  or  inconsistent  data 

entries  are  sensed  and  the  question  reissued  to  the  user.  The 

program  considers  up  to  ten  choices  of  molecule,  isotope,  vibrational 
band,  and  intensity  level  in  a  specified  frequency  range.  From  The 
selection,  output  can  be  generated  in  several  specific  forms. 

The  new  file  structure  of  HITRAN  is  shown  in  Figure  !.  A 
complete  description  of  the  changes  and  improvements  to  HITF  » v  a 
well  as  sources  of  the  data,  band  sums,  etc.  can  be  found  " 

3.  Table  6  illustrates  the  actual  image  of  the  new  dai. .  ■  ■  ;oi 

individual  quantum  transitions  for  a  small  frequet.v  '  <ing.;  ;  h 
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format  c:.h.f::.6.if:eio 

-  100  characters  pe 


3  . OP2F5 . - , F10  a  .  F- . 2 , F8 . 5 
transition 


3.2A9. 311. 312' 


thi:  format  correspond  .  as  follows: 


Mo 


iso 


*■- 

c 


7 


n 

v 


12  - 


E1C.3 

E1C.3 

F5  .  - 

F5 .  - 

F10.4 

F4.2 

F8.5 


IER  coce 


0 

1 


u 


6 


Molecule  number 

Isotope  number  (1-  most  abundant.  2-  second,  etc.) 

Freaue-  tv  in  cm 

Inters. tv  in  cm"  /(molec-cm  )  __ 2 9 6 K 

-  Transition  probability  in  Debves'  (presently  lacking  internal  partition,  sht 
Air  -  broadened  halfwidth  (HVHM)  in  cm  |atm  (2  2  9  6  K 
Self-broadened  halfwidth  ^HwllM)  in  cm"  /atm  g  296H 
Lower  state  energy  in  cm 

Coefficient  of  temperature  dependence  of  air-broadened  halfvirith 

Shift  of  transition  due  to  pressure  (presently  empty;  some  coupling  coeffic:*: 

inserted) 

Upper  state  global  quanta  index 
Lower  state  global  quanta  index 
Upper  state  local  quanta 
Lower  state  local  quanta  , 

Accuracy  indices  for  frequency  ,  intensity,  and  halfwidth 

Indices  for  lookup  of  references  for  frequency,  intensity,  and  halfwi.itn  ( net 
presently  used) 

for  frequency  when  used: 

estimated  error  in  wavenumber 
>  1 . 

>  0.1  anc  <  1.0 

>  0.01  and  <0.1 

>  0.0C1  and  <  0.01 

>  0 . 0001  and  <  0 . 001 

>  0.00001  and  <  0.0001 
<  0.00001 


Table  6.  Example  of  Direct  Xma^e  of  Lines  on  tne  1986 
HITRAII  Database 
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first  file  on  the  distribution  magnetic  tape  is  the  program  SELECT 

discussed  above.  The  file  is  in  a  "card"  image  form,  i.e.  80  characters 
per  record  and  is  followed  by  an  end-of-file  marker.  The  program  is 
written  in  FORTRAN  77,  follows  all  ANSI  standards  and  should  run  on 
all  systems  that  have  FORTRAN  compilers. 

The  second  file  on  the  distribution  tape  is  the  database. 
There  are  5100  characters  per  record  for  this  file  followed  by  an 

end-of-file  marker. 

The  third  file  on  the  tape  is  presently  a  table  of  the 
isotopic  abundances  of  the  species  on  the  database,  a  copy  of  this 
table  is  given  in  Table  2.  In  the  future,  this  file  will  also  contain 
reference  tables,  partition  sums  and  information  of  this  nature.  This 

file  contains  80  characters  per  record  and  ends  with  an  end-of-file 

marker. 


The  fourth  and  last  file  on  the  tape  is  for  heavy  trace 
constituents  of  the  atmosphere  that  have  dense  spectra  comprised  of 
many  thousands  of  lines.  In  the  last  file  are  absorption  cross- 
sections  for  these  species  to  be  used  in  simulation  codes  such  as 
FASCODE.21  The  data  is  from  the  experimental  work  of  Massie  et 
al.22  The  structure  of  the  file  is  as  follows,  a  100  character  header 
card  containing  the  species  name,  the  number  of  cross-sections,  the 
start  and  end  frequency,  and  the  experimental  conditions  used  to 
record  the  data  comes  first.  Next  the  absorption  cross-sections  are 
given  in  10E10.3  format  per  line  and  the  first  entry  corresponds  to 
the  start  frequency  and  for  each  following  entry  one  must  add  the 
frequency  interval  defined  by 

_  (end  frequency  -  start  frequency) 
number  of  points  for  the  species 

This  is  repeated  for  each  species  given  in  Table  7.  The  file  is  blocked 
5100  characters  per  physical  tape  record  (same  as  the  database) 
with  the  end-of-file  marker  at  the  end. 
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Table  7.  Approximate  Absorption  Cross-Sections  on  File  4  of  HITRAN  Tape 


Species 

#  Cross- 
Sections 

Frequency 

Minimum 

cm' 1 

Maximum 

CIONCb 

5020 

765.0020 

819.9980 

HNO4 

5476 

770.0070 

829.9990 

F21 

5020 

7  8  5.0000 

839.9950 

CCI4 

1826 

786.0010 

805.9980 

FI  1 

2738 

830.0090 

859.9990 

F 1 2 

7301 

860.0080 

939.9960 

FI  14 

1  2320 

1025.0090 

1 159.9920 

F21 

4563 

1050.0040 

1099.9920 

FI  1 

3651 

1060.0040 

1099.9980 

F 1  2 

10039 

1070.0050 

1 179.9940 

FI  13 

12777 

1090.0080 

1229.9980 

FI  14 

11771 

1 160.0250 

1288.9920 

N2O5 

4647 

1225.0010 

1265.0000 

HNO3 

7301 

1270.0040 

1349.9930 

CIONO2 

3650 

1270.007 

1309.9910 

CF4 

1095 

1275.0030 

1286.9900 

N2O5 

11616 

1680.0030 

1780.0000 
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To  aid  us  in  keeping  track  of  all  the  updates,  additions 
and  changes  that  were  to  be  done,  an  agenda  system  was  created. 
The  system  was  based  on  the  dBASE  III  database  language  and  set 
up  for  an  1BM-PC  environment.  The  agenda  system  allowed  us  to 
list,  track,  modify  and  add  update  information.  Each  record  pertains 
to  a  proposed  update  from  a  single  source.  The  structure  of  the 
records  is  given  in  Table  8. 

Several  dBASEIII  programs  have  been  written  to  make 
the  database  Agenda  easier  to  use.  The  first  is  a  custom  menu 
editing  program.  When  the  command  "DO  CUMENU"  is  entered  in 
dBASEIII  the  screen  shown  in  Figure  2  appears.  To  modify  existing 
records  a  2  is  entered  and  the  screen  shown  in  Figure  3  is  displayed. 
The  fields  of  the  Agenda  are  shown  in  reverse  video.  Changes  are 
made  by  typing  in  the  corrections,  moving  through  the  various  fields 
is  accomplished  by  the  numeric  keypad  arrows.  Addition  of  a  new 
record  (enter  1)  gives  a  screen  similar  to  Figure  3  with  no 
information  in  the  reverse  video  area.  The  information  is  entered, 
followed  by  a  carriage  return  or  when  the  field  is  full  the  cursor 
automatically  moves  to  the  next  field  to  be  entered.  The  use  of  a 
custom  screen  greatly  simplifies  the  maintenance  of  the  agenda. 

The  other  feature  we  have  added  for  the  agenda  is  a 
print  program.  By  entering  the  command  "DO  LISTMOL"  the 
following  prompt  is  seen: 

ENTER  WHICH  MOLECULE  TO  LIST 

TYPE  ALL’  TO  GET  ENTIRE  LIST 

MOLECULE  = 

by  typing  ALL  the  entire  data  base  is  listed  three  records  per  page. 
To  get  all  records  pertaining  to  a  particular  molecule  enter  the 
chemical  symbol,  e.g.  for  water  enter  H20,  peroxide  H202,  etc.  A 
second  program,  LISTMOL2,  creates  a  file  called  AGENDA. LST  which 


Number  of  data  records:  80 
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LOCATION  Character  80  Location  of  data  file 

FILENAME  Character  40  Name  of  permanent  file 

STATUS  Character  240  Description  of  the  state  the  data  is  in 

COMPLETION  Date  8  Date  data  is  added  to  atlas 


EXIT  AGENDA 


DO  YOU  WISH  TO 

1 .  ADD  A  RECORD 

2.  MODIFY  EXISTING  RECORD 

3.  QUIT 

ENTER  YOUR  CHOICE  AND  PRESS  ENTER 
Press  any  key  to  continue  .  .  . 

Figure  2.  Custom  Editor  Menu 


ENTER  MOLECULE  C02 
ENTER  ISOTOPE  111  11111 


COMPLETION  11/30/85 


FREQUENCY 

HW-AIR 

ENERGY 

SHIFT 

VLOWER 

VIBRATION 

LOCATION 

FILE  NAME 

SOURCE 


y  STRENGTH 

y  HW-SELF 

y  TDEP-HW 

n  USE 

400  VUPPER 

All  bands 
File  on  VAX 
C02:C02L1NES.DAT 
Rothman,  Global  analysis 


y 

n 

n 

y 

95  00 


STATUS  Needs  to  have  03301-00001  data  from  Rinsland  added; 

needs  to  have  perturbed  bands  added.  Needs  to  have 
self-broadened  halfwidth  and  temp-dependence  added. 


TO  EXIT  PRESS  PgDn  KEY 


Figure  3.  Screen  for  Editing  Records 


can  be  transferred  to  the  Cyber  or  Vax  systems  to  be  printed  by  the 
laser  printers.  The  form  is  three  records  followed  by  a  page  eject. 

We  now  have  an  automated  means  of  keeping  track  of 
changes  (post,  present,  and  future)  made  to  the  database. 
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2.0  TEMPERATURE  DEPENDENCE  OF  THE  HALFWIDTH 


It  became  apparent  during  the  literature  search  that  the 
temperature  exponents  for  the  halfwidths  are  not  well  understood. 
The  values  used  on  the  GEISA  atlas  and  many  of  the  HITRAN  values 
come  from  studies  of  only  a  few  transitions.  We  undertook  the  task 
of  evaluating  the  temperature  dependence  of  the  halfwidths  for  two 
asymmetric  rotors,  H2O  and  O3.  In  our  studies  we  looked  at  both  the 

vibrational  and  rotational  structure  of  this  parameter. 


Temperature  Dependence  of  Ozone  Halfwidth; 


The  temperature  dependence  of  N2-broadened 
halfwidths  of  ozone  was  studied  over  a  temperature  range  of  200- 
500K.  We  assumed  the  optical  cross-section  to  vary  as  a 
temperature  independent  part  times  temperature  to  a  power  m  (o(T) 
=  00  •  Tm).  This  relates  the  halfwidth  at  temperature  T  to  the  value 
at  a  reference  temperature  Tq  as 


Y(T)  =  *T0) 


ll°\" 

\T/ 


(10) 


The  exponent  n  was  determined  by  calculating  the  halfwidth  for  126 
ro-vibrational  transitions  of  ozone  at  seven  temperatures.  The 
halfwidths  were  calculated  using  quantum-Fourier-transform  theory 
with  improved  dynamics  (QFT-ID).2^  The  transitions  chosen 
considered  a  wide  range  of  rotational  quantum  number,  J"  (1  to  35), 
the  full  manifold  of  Ka”  for  a  particular  J",  and  transitions  important 
to  stratospheric  studies.  The  vibrational  states  considered  both  A- 
and  B-type  transitions  using  constants  for  the  ground  state  and  V3 
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bands.  All  molecular  constants  were  considered  the  best  available  at 
the  time  of  the  calculations. 

For  each  transition  studied  linear  regression  was  applied 
to  the  coordinates;  ln(T2/T  j )  and  ln(y(T \  )/y(T2)),  where  T2,  the 
reference  temperature  was  taken  at  296K.  From  linear  regression 
we  obtain  the  slope  n,  i.e.  the  temperature  exponent  and  the 
correlation  coefficient  of  the  fit.  Figure  4  shows  such  a  plot  for  the 
27  15  12  <-  27  13  15  transition  of  ozone.  For  all  the  transitions 
studied,  the  correlation  coefficients  were  equal  to  one  to  four 

significant  figures  comfirming  the  assumed  temperature  of  the 
optical  cross-section,  i.e.  oo#Tm. 

The  results  of  the  calculations  give  an  average 

temperature  exponent  of  n  =  0.77  for  nitrogen-broadening.  The 

temperature  exponents  are  transition  dependent  but  the  fluctuation 
about  the  average  value  is  small  (±6%).  The  value  of  n  as  a  function 
of  J",  and  Ka"  was  investigated  and  the  work  was  related  to  other 
work  in  the  literature.  A  complete  description  can  be  found  in 

reference  15  or  Appendix  B.  This  work  lead  to  an  average 
temperature  exponent  n  =  0.76  for  air-broadening  of  ozone. 


2.2  Temperature  Dependence  of  Halfwidths  of  Water  Vapor 

We  have  performed  calculations  to  study  the 
temperature  dependence  of  ^-broadened  halfwidths  of  H2O.  In  the 
calculations  we  employed  cutoff-free  theory  with  improved 
dynamics  and  ATC  interruption  functions.  Some  50  water  vapor 
transitions  were  considered  for  the  pure  rotation  and  V2  bands.  The 
transitions  were  supplied  to  us  by  Dr.  Wesley  Traub  of  the  Harvard- 
Smithsonian  Center  for  Astrophysics  and  represent  some 
stratospherically  important  lines  of  water  vapor.  Nine  temperatures 
were  considered  from  200  to  1000  K  to  allow  the  results  to  be 
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applicable  to  the  stratosphere  as  well  as  hot  sources.  The 
calculations  were  done  and  the  results  are  given  in  Table  9.  One 

surprising  fact  was  the  similar  results  for  the  pure  rotation  and  V2 
bands.  We  had  expected  more  of  a  dependence  on  the  vibrational 

band.  In  general,  the  V2  results  for  n  were  lower  than  the  pure 

rotation  n  values.  However,  there  were  some  exceptions  to  this.  The 
large  variation  in  n  with  the  rotational  quantum  numbers  was  also 
surprising.  Average  values  were  calculated  for  each  band  giving  n 
(pure  rotation)  =  0.684  with  a  +15  to  -28%  variation  about  the 
average  and  n  (V2)  =  0.676  with  a  +24  to  -10%  fluctuation  about  the 
average.  This  can  be  compared  with  the  value  we  calculated  for  the 
ozone-nitrogen  system  of  n  =  0.77  with  a  ±6%  fluctuation.  It  is 
suggested  that  a  value  of  n  =  0.68  be  used  for  water  vapor  and, 

whenever  possible,  use  a  calculated  value  from  Table  9.  A  complete 
description  of  the  work  is  given  in  reference  24  or  Appendix  B. 


to  to  to 


Table  9.  Temperature  Exponents  for  H2O-N2  System 


Assuming  the  form. 


W=y(t„) 


/Mn 
\t  I 


Transition 

J"  Ka"  Kc"  J'  Ka*  Kc' 


Pure  Rotation  V2 

n  n 


1112  0  2 

1  10  2  2  1 

113  2  2 

2  13  3  0 

2  0  3  3  1 

3  0  3  3  3  0 

3  1  3  4  2  2 

3  12  4  2  3 

3  2  14  3  2 

3  3  1  4  4  0 

3  3  0  4  4  1 

4  0  4  4  3  1 

4  0  4  5  15 

4  14  5  0  5 

4  1  4  5  2  3 

4  13  4  4  0 

4  2  3  5  14 

4  2  3  5  3  2 

4  2  2  5  3  3 

4  3  2  4  4  1 

4  3  1  4  4  0 

5  0  5  6  1  6 

5  0  5  5  3  2 

5  15  4  0  4 

5  15  6  0  6 

5  14  5  4  1 

5  1  4  5  2  5 

5  2  4  6  15 

5  2  3  6  3  4 

5  3  3  5  4  2 


.7  1  82 
.7259 
.7086 
.6786 
.6713 
.7108 
.7394 
.7067 
.6964 
.6507 
.6585 
.7133 
.6859 
.6729 
.7406 
.7045 
.7064 
.7169 
.7171 
.6696 
.6937 
.6106 
.7085 
.6859 
.5939 
.7037 
.7036 
.6682 
.7274 
.6765 


.7110 

.7067 

.7032 

.6708 

.6737 

.7130 

.7441 

.7074 

.6988 

.6577 

.6641 

.7084 

.6841 

.6664 

.7435 

.7134 

.6997 

.7003 

.7212 

.6643 

.6969 

.6076 

.6954 

.6840 

.5871 

.7052 

.7060 

.6590 

.7301 

.6581 


5  3  3  6  2  4 

5  3  2  5  4  1 

6  0  6  6  15 

6  1  6  5  4  1 

6  16  6  2  5 

6  2  5  5  5  0 

6  2  5  6  3  4 

6  3  4  6  4  3 

6  3  3  6  4  2 

6  5  2  6  6  1 

6  5  16  0  6 

8  5  3  8  5  4 

9  3  6  9  5  4 

9  4  5  8  5  4 

10  5  5  1  0  4  6 

11  4  8  11  5  7 

12  6  6  12  4  8 

14  9  514  8  6 


.7277  .7284 

.7258  .7  169 

.6407  .6389 

.6262  .6124 

.5833  .5689 

.6472  .6370 

.6378  .6120 

.6619  .6430 

.7498  .7344 

.7192  .7348 

.6209  .6211 

.6477  .6580 

.7087  .6979 

.7228  .7058 

.7  1  1  3  .6858 

.5868  .5860 

.6441 

.5170 
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3.0  WATER  VAPOR  DATA  IN  THE  VISIBLE  SPECTRAL  REGION 


We  have  taken  the  data  for  water  vapor  in  the  visible 
region  of  the  spectrum  from  Camy-Peyret  et  ai. 25-26  The  data 
consist  of  4612  lines  in  the  region  13,435  cm-'  to  22,597  cm'1.  This 
compares  with  2932  water  vapor  lines  presently  on  the  database  in 
this  range.  The  data  contains  the  bands  listed  in  Table  10  and 
includes  eleven  bands  that  were  not  present  on  the  1982  atlas. 

The  intensities  are  reported  in  units  of  cm'2  atm'1  at  300 
K  and  were  converted  to  the  units  used  on  the  database,  cm'l/molec 
cm *2  and  to  296  K,  the  reference  temperature  of  the  database.  This 
first  is  accomplished  by  dividing  by  Loschmidt's  constant,  Nj  = 
(Na/Vm)  =  2.68675  x  lO1^  cm'3«moIec,  corrected  to  300  K  by  the 
Ideal  Gas  Law,  i.e. 

NL(300K)  =  NL(273.15K)*(273. 15/300). 

After  which  the  temperature  correction  must  be  done. 

The  data  file  contains  the  wavenumber  of  the  transition 
and  the  corresponding  error,  the  measured  intensity  and  estimated 
uncertainty,  the  vibrational  and  rotational  quantum  numbers,  and 
the  absorption  at  the  line  center  in  tenths  of  a  percent.  Many  of  the 
lines  in  the  data  do  not  contain  the  intensity  and  some  are 
unassigned  lines.  The  data  file  was  split  into  lines  with  the  measured 
intensity,  1319  lines,  and  the  lines  without  a  measured  intensity, 
3293  lines.  The  lines  that  have  intensities  can  be  put  into  the  new 
format  without  difficulty.  Those  without  intensities  cannot  be  used 
without  some  estimation  of  the  intensities  for  the  lines.  Fortunately 
the  absorption  at  the  line  center  is  given  and  can  be  related  to  the 
intensity  of  the  line. 


Table  10.  Vibrational  Bands  Present  in  the  Water  Vapor  Data  from 


Camy-Peyret  et  al. 


Vibrational  Bands 


Bands  present  in 

1982  Atlas 

221, 

301, 

202, 

042, 

320, 

141,  023 

Bands  New  to  the 

331, 

213, 

312, 

411, 

501, 

313, 

Database 

303, 

403, 

402, 

421, 

412 

For  the  lines  with  intensities,  the  intensity  was  converted 
to  the  AFGL  units  by  division  by  Loschmidt's  number  at  300  K.  The 
temperature  correction  of  the  intensity  to  296  K  was  then  done.  The 
file  was  reformatted  into  the  1986  format  and  the  halfwidths  of 
Gamache  and  Davies27  added  to  the  lines.  The  error  estimates  for  the 
data  were  added  to  the  new  iormai  in  the  IER  field.  The  values  used 
for  IER  are  listed  in  Table  11. 

For  lines  that  only  have  peak  absorption  values  the 
intensities  must  be  estimated  by  a  model.  A  simple  Lorentz  model  of 
the  line  shape  was  used  to  relate  the  peak  absorption  values  to  the 
line  intensity. 

Of  the  lines  of  data  that  have  no  measured  intensities, 
1842  of  these  have  quantum  assignments  (this  leaves  1451 
unassigned  lines).  For  all  the  lines  the  peak  absorption  at  the  center 
of  the  line  is  given.  The  intensities  for  these  lines  can  be  estimated 
from  the  peak  absorption  values  for  the  lines  that  have  measured 
intensities  and  a  model  that  relates  peak  absorption  to  the  measured 
intensity.  The  model  of  line  absorption  chosen  for  this  was  the 
Lorentz  line  shape.  The  absorption  coefficient  for  the  Lorentz  line 
shape  is  given  by 


Where  S  is  the  integrated  line  intensity  and 
Y  is  the  Lorentz  halfwidth  and 
vo  is  the  frequency  at  line  center. 

The  Lorentz  shape  has  the  following  form. 


Table  11.  1ER  Values  for  Wavenumber  and  Intensity 


IER  Value 


Estimated  Error  in  Wavenumber 


0 

1 

2 

3 

4 

5 

6 


>  .99999 

>  .  1  and  <  .9 

>  .01  and  <  .1 

>  .001  and  <  .01 

>  .0001  and  <  .001 

>  .00001  and  <  .0001 
<  .00001 


IER  Value 


Estimated  Error  in  Intensity 


0 

1 

2 

3 

4 

5 

6 
7 


<  l.% 

1  %  anti  <  2.5%, 

>  2.5 %  and  <  5.0% 

>  5.0%  and  <  10.0% 

>  10.0%  and  <  15.0%, 

>  15.0%  and  <  20.0% 

>  20.0% 

Estimated  Intensities 
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Peak  absorption  at  the  center  of  the  line  can  be 
associated  with  S/rcy.  From  the  1319  lines  of  data  that  contain  both  S 
and  peak  absorption,  P,  at  the  line  center,  we  can  equate  the  two  via 
a  scale  factor,  i.e. 


(12) 


In  the  calculations,  the  factors  of  n  and  y  have  been  absorbed  into  the 
scale  factor  so  for  each  of  the  1319  lines  f  =  S/P  is  calculated  for  each 
line.  Before  the  scaling  was  done  the  intensities  were  converted 
from  cm~2  atm*1  units  to  the  units  used  by  AFGL,  cm1/  molec  cm*2. 
The  average  value  of  the  scale  factor  was  calculated  from  the  data 
giving  f  =  0.583  x  10*25,  Multiplying  the  peak  absorption  of  a  line  by 
f  will  give  the  estimated  strength  at  300  K  in  AFGL  units.  This  was 
tested  by  recalculating  the  intensity  for  the  1319  lines  using  f  and 
comparing  these  values  with  the  measured  strengths. 

This  approach  did  not  work;  the  average  error  from  the 
method  was  approximately  2000.  The  reason  for  this  was  that  three 
lines  had  a  scale  factor  of  order  10*22  to  10*23,  whereas  all  the  other 
lines  were  of  order  10*26,  Hence  the  average  factor  essentially 
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came 


from  the  three  large  factors.  To  correct  for  this  the  three 
troublesome  lines  were  removed.  Then  the  standard  deviation  of  the 
scale  factor  was  calculated  and  then  a  scale  factor  was  calculated 
using  only  factors  within  ±1  a  of  the  first  average.  This  gave  (note 
three  troublesome  factors  were  removed) 

#  #  Lines  ave  f 

Lines  ave  f  a  within  1  a  within  1  a 

1316  0.408  x  10'26  0.315  x  10*26  1259  0.356  x  10'26 

The  new  value,  f  =  0.356  x  10*2^  was  then  used  to 
recalculate  the  strengths  and  compared  with  the  experimental  values 
for  the  1319  measured  lines  (note  three  "troublesome"  lines  included 
here)  the  results  were 

Average  %  Dif  =  -  1.353 

Average  Absolute  %  Dif  =  16.912 

These  results  suggest  the  model  can  supply  the  line 
intensities  for  the  lines  where  only  a  peak  value  is  given.  The  lines 
of  data  without  a  line  intensity  were  taken  and  the  intensity 
estimated  via  the  model,  i.e.  S  =  fP  with  f  =  0.356  x  10'2^.  This 
produced  a  line  file  (at  300  K)  that  can  be  incorporated  into  the  AFGL 
database  after  rescaling  to  296  K. 


The  simple  Lorentz  line  shape  model  was  discussed  with 
Jean-Marie  Flaud,  University  of  Paris.  He  has  suggested  that  a 
different  model,  based  on  the  curve-of-growth  (COG)  method,  be 
tried.  In  the  COG  model,  the  peak  absorption  is  given  by 


P  =  ^(o0)=  1  -e(  k°N  r0^Ou) 


which  assumes  a  Voigt  profile, 
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1 

log  2 

Yd 

71 

A 


2  -  -  lAl t 


k(0,  log  2) 


Vd 


(14) 


If  one  assumes  that  the  variation  of  <> v ( 0 0 )  *s  mainly  due  to  the 
variation  of  Yd  which  is  proportional  to  oq  one  has 


0jco)  =  — ,  a  =  effective  constant 

Thus  equation  13  can  be  written  as 
P  =  1  -e(k>a/^ 

and  with  the  definition  (3  =  Nla  the  following  can  be  obtained 


(15) 


(16) 


N  C7r 

ka= - -  ln(l-P). 

P  (17) 

From  the  values  of  line  intensity  and  peak  absorption,  the  constant  P 
is  evaluated  after  which  the  intensities  can  be  estimated  for  lines 
that  have  only  peak  absorption  values  using  the  average  P  value.  In 

the  data  set  there  are  1311  lines  that  have  intensities  and  peak 

absorption  values  reported.  The  effective  constant  p  has  been 
evaluated  for  these  lines.  At  the  suggestion  of  Flaud,  we  have 

averaged  p  for  the  two  regions  13500-16500  cm'1  and  16500-23000 
cm'1  (see  Table  12).  An  average  over  all  lines  was  done  and  an 
average  as  a  function  of  P  for  the  ranges  shown  in  Table  12.  For 
peak  absorption  values  greater  than  150  the  average  p  values  are 
similar  for  both  regions.  For  the  values  less  than  150  there  is 
considerable  change  for  the  lines  greater  than  and  less  than  16500 

cm'1,  especially  for  the  peak  absorption  lines  between  75  and  150  in 
the  upper  wavenumber  region. 
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Table  12.  Average  Values  of  the  Effective  Constant 
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The  data  indicate  a  different  average  for  the  two  regions 
and  this  was  used  in  the  generation  of  the  intensities  for  the  lines 

with  peak  absorption  only.  Before  estimating  these  intensities,  the 
intensities  of  the  data  with  measured  intensities  was  calculated  via 
equation  17  with  the  average  values  of  [5.  [3  =  0.100715E9  for  region 
1  and  (3  =  0.22147E9  for  region  2.  The  average  absolute  percent 
error,  AAPD,  the  average  error,  E,  and  the  maximum  error,  Emax 
were  computed  using  the  measured  intensities  for  the  1311  lines. 
The  results  are  the  AAPD  =  13.0%,  E  =  -0.00073%  and  Emax  =  845%, 
this  is  compared  with  ihe  results  of  the  Lorentz  model,  AAPD  =  17% 
and  E  =  -1.4%  (Emax  was  not  calculated  in  this  model).  It  appears 
that  the  curve  of  growth  model  gives  better  agreement  with  the 

experimental  intensities  than  the  Lorentz  model  does.  Thus  the 
visible  spectrum  of  water  vapor  was  recalculated  via  the  COG  model 
in  the  same  manner  that  was  done  for  the  Lorentz  model.  Using  the 
COG  model  the  intensities  were  calculated  for  lines  with  only  the 
peak  absorption  given.  These  lines  were  merged  with  the  data 
containing  both  S  and  P  values.  The  data,  which  is  for  300  K,  was 
then  corrected  for  296  K,  the  reference  temperature  of  the  database. 
A  similar  procedure  was  used  on  the  data  without  quantum 
identifications  (unassigned  lines).  This  yields  a  file  with  intensities 
for  all  the  lines  reported  by  Camy-Peyret  et  al.25-26 

The  water  vapor  data  for  the  visible  region  of  the 
spectrum  is  almost  complete.  The  final  step  is  to  correct  the 
intensities  of  the  saturated  transitions,  i.e.  those  that  had  peak 

absorptions  greater  than  500.  Jean-Marie  Flaud  has  sent  a  list  of  the 

improved  intensities  for  these  lines,  see  Table  13.  To  make  the 
correction,  the  transitions  were  selected  from  the  original  water 
vapor  data  in  Flaud's  units  of  cm'^atnv1.  The  intensities  were  then 
corrected  to  the  values  in  Table  13  and  the  file  was  then  converted 
to  AFGL  units  and  format  and  then  corrected  from  300  K  to  296  K. 
The  final  correction  file  is  given  in  Table  14.  The  file  was  merged 
with  the  other  lines  for  the  visible  region. 
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Table  13.  Flaud's  Notes  for  Lines  with  P  >  500 


cm "2  atm'l  at  300  K 


13628.8617 

000 

-  101 

22  1 

552 

-  4.74-5 

13628.91  13 

542 

-  643 

301 

488 

-  3.52-5 

13631 .9693 

717 

-  818 

301 

866 

-  1.12-4 

1  3690.0014 

542 

-  541 

22  1 

511 

-  3.63-5 

13695.0327 

440 

-  441 

22  1 

802 

-  8.60-5 

13708.5954 

551 

-  550 

301 

857 

-  1.05-4 

404 

-  515 

202 

13766.4625 

551 

-  550 

22  1 

792 

-  8.09-5 

13768.9036 

524 

-  523 

301 

642 

-  5.56-5 

13768.9739 

642 

-  643 

301 

551 

-  3.79-5 

13782.1900 

441 

-  440 

301 

800 

-  8.39-5 

13782.2523 

440 

-  441 

301 

995 

-  2.44-4 

13893.1820 

432 

-  331 

301 

627 

-  5.30-5 

13897.4973 

312 

-  211 

301 

992 

-  2.54-4 

13914.8369 

704 

-  6.13-5 

13914.9602 

413 

-  312 

301 

655 

-  6.13-4 

13915.0509 

524 

-  423 

301 

999 

-  3.48-4 

13929.3493 

707 

-  606 

301 

889 

-  1.23-4 

13942.5475 

909 

-  808 

301 

512 

-  3.90-5 

845 

-  744 

301 

13947.2608 

624 

-  523 

301 

983 

-  2.19-4 

15205.2919 

515 

-  616 

3  1  1 

531 

-  3.96-5 

15402. 1988 

313 

-  212 

3  1  1 

592 

-  4.97-5 

15437.8123 

413 

-  312 

3  1  1 

558 

-  4.46-5 

16888.2383 

220 

-  221 

203 

563 

-  6.50-5 

B  =  blended 
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P  =  perturbed 


0-0-  CQ  Cu  0-  ft-  CQ 


13628 . 861700  0 . 198E-23 
1  3628.  911  300  0  .  154E-23 
13631.949300  0.490E-23 
1  3690  .  001  400  0 . 1S7E-23 
13695 . 032700  0 . 371E-23 
13708 . 595400  0 . 459E-23 
13  708.59  5400  0. 4  4  7E-2  3 
13766 .462500  0 . 354E-23 
1  3  768  .  903600  0 . 238E-23 
13768.973900  0.166E-23 
1  3782  .  190000  0  .  36  1E-23 
13782.252300  0. 105E-22 
13893.182000  0.226E-23 
13897 . 497300  0 . 107E-22 
13914.836900  0.256E-24 
13914.960200  0.259E-22 
13915.050900  0.148E-22 
13929.349300  0.528E-23 
13942.547500  0.170E-23 
1  394  2  .  547500  0  .  172E-23 
13947.260800  0.940E-23 
15205.291900  0.170E-23 
15402.198800  0.208E-23 
15437.812300  0.188E-23 
16888  .  238  300  0  .  274E-23 


2379400. 0000 

0 

0 

0 

1 

0 

1 

2 

2 

1 

c 

0 

0 

16  . 

756.72500.0000 

5 

•3 

2 

6 

4 

3 

3 

0 

1 

0 

0 

0 

161 

744 .  1  6  300  .0000 

7 

i 

7 

8 

1 

8 

3 

0 

1 

0 

0 

0 

161 

610 . 34100 .0000 

5 

4 

2 

5 

4 

1 

2 

2 

1 

0 

0 

0 

16  1 

488.10800 .0000 

4 

4 

0 

4 

4 

1 

2 

2 

1 

0 

0 

0 

161 

742 .07600 .0000 

5 

5 

1 

5 

5 

0 

3 

0 

1 

0 

0 

0 

16  1 

326 . 62500 .0000 

4 

0 

4 

5 

1 

5 

2 

0 

2 

0 

0 

0 

161 

742 . 07600  .0000 

5 

5 

1 

5 

5 

0 

2 

2 

1 

0 

0 

0 

161 

446 . 51100 . 0000 

5 

2 

4 

5 

2 

3 

3 

0 

1 

0 

0 

0 

161 

756  . 72500  .0000 

6 

4 

2 

6 

4 

3 

3 

0 

1 

0 

0 

0 

161 

488.13400  .0000 

4 

4 

1 

~i 

4 

0 

3 

0 

1 

0 

0 

0 

161 

488 . 1  3  4  0  Q  .0000 

4 

4 

0 

4 

4 

1 

3 

0 

1 

0 

0 

0 

161 

28  3 . 21900 .0000 

4 

3 

2 

3 

3 

1 

3 

0 

1 

0 

0 

0 

161 

95.17600. 0000 

3 

1 

2 

2 

1 

1 

3 

0 

1 

0 

0 

0 

161 

0 .00000 .0000 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

161 

1  73 . 36500 .0000 

4 

1 

3 

3 

1 

2 

3 

0 

1 

0 

0 

0 

161 

300 . 36200 .0000 

5 

2 

4 

4 

2 

3 

3 

0 

1 

0 

0 

0 

161 

446 .69700 .0000 

7 

0 

7 

6 

0 

6 

3 

0 

1 

0 

0 

0 

161 

744.06400.0000 

9 

0 

9 

8 

0 

8 

3 

0 

1 

0 

0 

0 

161 

927.74400  .0000 

8 

4 

5 

7 

4 

4 

3 

0 

1 

0 

0 

0 

161 

446  .  SHOO  .0000 

6 

2 

4 

5 

2 

3 

3 

0 

1 

0 

0 

0 

161 

447 .25200 .0000 

5 

1 

5 

6 

1 

6 

3 

1 

1 

0 

0 

0 

161 

79  .  49600 .0000 

3 

1 

3 

2 

1 

2 

3 

1 

1 

0 

0 

0 

161 

173.  36500. 00 '.'0 

4 

1 

3 

3 

1 

2 

3 

1 

1 

0 

0 

0 

161 

1  34.  90200  .0000 

2 

2 

0 

2 

2 

1 

2 

0 

3 

0 

0 

0 

1  6  1 

Table  iu.  Corrected  II0C  Lines  fcr 


which  F  >  500 


Before  adding  the  data  to  the  database  and  publishing  the 
results  we  are  going  to  study  the  data  (e.g.  FASCODE  plots,  etc.)  and 
compare  the  results  with  other  data  for  this  region. 
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4.0  MISCELLANEOUS  TASKS 


Throughout  the  course  of  the  contract  there  were  many 
tasks  that  were  performed  in  conjunction  with  other  groups  as  well 
as  preparing  data  for  the  next  version  of  the  database.  These  tasks 
are  discussed  below. 


4.1  Ozone  Data  for  Minor  Isotonic  Species 

A  file  of  lines  for  the  less  abundant  isotopic  species  of 
ozone  was  put  on  the  VAX  computer.  This  file  is  to  be  added  to  the 
database  to  upgrade  the  information  for  these  isotopes  of  ozone.  Not 
present  in  the  file  (1982  parameters)  is  the  halfwidths  for  the 
transitions.  The  file  consists  of  11235  lines  of  data  from  934  cm1  to 
1167  cm"'.  The  file  was  sent  to  the  CYBER  where  halfwidths  of 
Gamache  and  Rothman28  were  added  using  the  CYBER's  mass-storage 
facilities.  The  file  was  transferred  back  to  the  VAX  where  it  was 
reformatted  to  the  1986  format  and  the  additional  parameters  (e.g. 
n,  R2,  etc.)  added.  The  file  is  now  on  the  VAX  as  A86:OZONE_ISO.DAT. 


4.2  Ozone  Halfwidths  for  the  GEISA  Atlas 

A  magnetic  tape  containing  the  ozone  data  that  is  present 
on  the  GEISA  atlas  was  received  from  Nicole  Husson.  The  tape 
contained  two  files,  one  containing  the  temperature  exponent  for  the 
halfwidth  and  the  other  file  has  a  GEISA  date  code  in  this  field. 
Because  we  have  the  software  readily  available  to  add  our  ozone 
halfwidths28  to  the  GEISA  atlas,  we  offered  to  do  this  task  for  Dr. 
Husson.  This  will  ensure  that  no  mis-communicalion  can  occur  and 
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that  the  GE1SA  halfwidths  for  ozone  are  correct  and  match  the 
HITRAN  results  line-for-line,  default  halfwidths  included. 


The  files  were  taken  from  the  tape  using  the  AFGL 
copany  routine.  The  tape  was  labelled  as  follows 

LABEL,TAPE1,VSN=M2764,CV=EB,LK=KU,D=PE,P0=R,F=L 

and  copany  needed  a  file  card  to  describe  the  77  card  images  per 
physical  record  on  the  tape,  i.e. 

FILE, TAPE  1,FL=80,RB=77,RT=F,BT=K. 

By  running  copany  twice  we  were  able  to  extract  each  file  and  save 
them  as  direct  access  files.  The  ADDHW03  program  was  modified  to 
read  the  GEISA  data  and  using  the  mass  storage  file  of  O3  halfwidths, 
add  them  to  the  output  file.  The  output  files  were  checked  to  ensure 
the  program  ran  as  expected  and  the  two  files  were  then  transferred 
to  a  magnetic  tape  in  the  form  that  they  were  sent  to  us.  The  tape 
was  dumped  and  the  results  appeared  correct.  We  then  sent  the 
tape  to  Dr.  Husson  for  addition  to  the  next  version  of  the  GEISA  atlas. 


4,3  Halfwidth  and  Temperature  Dependence  of  Halfwidths  for 
Database 


At  suggestion  of  Linda  Brown  of  JPL,  we  generated  a 
table  of  the  range  of  the  halfwidths  and  temperature  exponents  of 
the  halfwidth  for  each  molecule  on  the  database.  To  do  this  a 
program  was  written  that  read  the  unblocked  database  and  stored 
the  maximum  and  minimum  values  of  y  and  n  for  each  molecule. 
After  the  program  was  run  the  output  was  taken  and  put  into 
tabular  form  (see  Table  5). 


Dr.  Mary  Ann  Smith  of  NASA  Langley  Research  Center 
requested  a  listing  of  ozone  shift  values  for  lines  in  the  region  1070 
cm'l  to  1140  cm-'.  Using  SELECT  all  ozone  lines  in  this  region  were 

removed  from  the  database.  From  this  file  only  lines  with  J  <_  35 

were  kept,  since  the  shift  calculations  stopped  at  J  =  35.  The 
calculated  shift  values  were  recovered  from  backup  tape  and  put  on 
the  CYBER  and  a  program  was  written  to  add  the  corresponding  shift 
to  the  lines  of  data  based  on  the  rotational  quantum  numbers  (no 
vibrational  dependence  was  considered).  This  only  produced  a  small 
number  of  shift  values  for  the  ozone  lines.  The  reason  was  that  the 
calculation  of  the  shift  is  specific  to  the  initial  and  final  state 
rotational  quantum  numbers.  The  listed  shift  values  mostly  had  the 
quantum  numbers  reversed  (f->  i,  i->  f)  from  the  lines  Smith 
requested.  The  calculations  of  the  shift,  however,  have  the  shift  of 
the  individual  energy  levels  listed  which  allows  the  calculation  of 

both  Avjf  and  Avfj.  This  is  described  below. 

If  one  calculates  the  shift  in  the  wavenumber  for  the 
transition  from  level  A  to  B,  the  shift  for  the  transition  B  to  A  can  be 
calculated  from  the  shifts  in  the  individual  levels  A  and  B.  This,  of 
course,  neglects  vibrational  dependence  of  the  shift  which  may  be 
important  for  many  systems.  To  see  this  consider  the  following 

diagram, 
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In  the  unperturbed  case,  the  wavenumber  of  the 

transition  is  given  by 

v(cnr')  =  EB(cm_1)  -  E^fcm'^).  (18) 

When  the  levels  are  pressure  shifted  the  wavenumber  is 
vs(cm“l )  =  (Eb+AB)  -  (E,\+AA) 

=  Eb~Ea  +  (AB-AA) 

=  v  +  (AB-AA).  (19) 

The  line  shift  is  given  by  Av  =  (AB-AA).  When  the  reverse  transition 
B-»A  occurs  the  shift  is  Av  -  (AA-AB).  Because  of  the  sign  of  the  shift 
of  the  individual  levels,  Av(A->B)  does  not  necessarily  equal  -A  v 
(B->A).  For  the  cases  when  AA  and  AB  have  the  same  sign  Av(A->B)  = 
-Av(B->A),  when  one  level  shifts  up  and  the  other  down,  Av(A-»B)  * 
Av(B-»A).  In  these  cases  the  shift  for  the  reverse  transitions  must  be 

calculated  from  the  shifts  of  the  energy  levels.  These  shifts  were 

calculated  in  the  QFT-ID  calculations  for  O3;  thus  the  reverse 
transitions  can  be  calculated. 

The  shifts  for  the  reverse  transitions  were  calculated  and 
put  into  a  file  with  the  rotational  quantum  numbers,  the  halfwidth 


and  shift.  The  shift  data  was  then  added  to  the  file  of  ozone  lines. 
This  supplied  data  for  all  the  lines  in  the  listing.  The  data  was  then 
sent  to  Dr.  Smith  for  comparison  with  experimentally  determined 
shifts  for  ozone  lines. 

The  same  experiment  has  yielded  ozone  halfwidths  and, 
in  a  preprint,  Smith  stated  that  for  v  j  lines  our  calculations  are  a 
fairly  constant  6%  lower  than  the  experimental  values.  For  the  few 
V3  lines  that  she  has  measured  the  calculations  average  about  3% 
lower  than  experiment.  The  reason  for  the  improvement  for  the  V3 
lines  is  that  the  calculations  were  done  for  A-  and  B-type  transitions 
and  the  vibrational  constants  used  were  from  the  V3  and  pure 
rotation  bands.  Hence  the  calculated  halfwidths  for  V3  transitions 
contain  the  vibrational  dependence  explicitly  whereas  the  vj 
halfwidths  were  calculated  using  ground  state  constants  only. 

Drs.  Ira  Nolt  and  Mian  Abbas  requested  ozone  halfwidths 
for  some  twenty  transitions.  The  halfwidth  data  is  to  be  used  for 
retrieval  experiments.  The  data  were  for  lines  in  the  40  cm-'  to  70 
cm-1  region  of  the  spewtrum.  The  data  is  listed  in  Table  15  and  was 
sent  to  Dr.  Abbas  at  Drexel  University. 


4,5  Halfwidths  for  Water  Vapor  for  HOT  Gas  Atlas 

Dr.  John  Selby  of  Grumman  Corp.  sent  a  magnetic  tape  of 
transitions  for  water  vapor  which  are  to  make  up  the  AFGL  HOT  gas 
atlas.  The  data  on  the  tape  are  for  the  reference  temperature  1000 
K.  Selby  requested  that  data  be  added  to  the  tape,  namely  the  air- 
broadened  halfwidths  and  the  temperature  dependence  of  the 
halfwidths.  The  programs  to  accomplish  this  task  reside  on  the 
CYBER  computer  and  make  use  of  mass  storage  files.  The  tape  was 
loaded  on  to  the  CYBER.  This  was  attempted  by  two  different 
methods;  first  the  CYBER  COPANY  routine  was  used  and  while  reading 


Table  15.  Ozone  Halfwidths  sent  to  Dr.  M.  Abbas 


Rotational 

Halfwidths 

Quantum 

Numbers 

cm'Vatm 

7 

7 

1 

6 

6 

0 

.0677 

11 

7 

5 

10 

6 

4 

.0691 

9 

9 

1 

8 

8 

0 

.0677 

19 

8 

12 

18 

7 

1 1 

.0685 

10 

6 

4 

9 

5 

5 

.0688 

12 

7 

5 

1 1 

6 

6 

.0690 

13 

10 

4 

12 

9 

3 

.0690 

13 

12 

2 

12 

1 1 

1 

.0669 

17 

6 

12 

16 

5 

1 1 

.0682 

17 

7 

1 1 

16 

6 

10 

.0684 

18 

6 

12 

17 

5 

13 

.0682 

18 

7 

11 

17 

6 

12 

.0683 

20 

7 

13 

19 

6 

14 

.0682 

20 

6 

14 

19 

5 

15 

.0682 

21 

6 

16 

20 

15 

.0683 

27 

5 

23 

26 

4 

22 

.0696 
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the  tape  a  parity  error  was  encountered.  Next  a  BUFFER  IN  program 
was  used  and  again  a  parity  error  was  encountered  while  the  tape 
was  being  read.  I  then  took  the  tape  to  the  AFGL  VAX  computer  and 
accomplished  copying  the  tape  into  a  file  without  any  problems. 

The  VAX  file  SELBY. H20,  18072  lines  of  data,  was 
transferred  to  the  CYBER  via  the  hyperchannel.  The  program  to  add 
the  halfwidths  along  with  the  five  halfwidth  files  that  are  input  to  it 
were  put  on  the  CYBER  from  backup  tapes.  The  halfwidths  were 
added  to  the  file  (note  the  file  is  already  in  the  1986  format)  and  the 
temperature  exponent  n  =  0.64  was  added  to  ail  the  lines.  The  file 
was  hyperchanneled  back  to  the  VAX  as  SELBYHW  H20.  The  file  was 
visually  inspected  along  with  the  output  from  the  CYBER  addition  of 
halfwidths  and  it  appears  to  be  correct.  The  file  was  written  to  a 
magnetic  tape  and  sent  back  to  Selby. 


4.6  COo  Data  from  L.  Rosenmann 

Laurence  Rosenmann  from  Ecole  Centrale  des  Arts  et 
Manufactures  sent  a  Macintosh  diskette  containing  halfwidths  and 

temperature  exponents  for  CCb  broadened  by  CO2,  H2O,  N2<  an(*  02- 

The  file  was  written  in  the  MacWrite  application.  The  file  was  loaded 
into  a  Macintosh  SE  and  the  MacWrite  file  was  converted  !o  text 

form.  lhe  tile  was  then  transferred  to  the  AFGL  VAX  computer  and 
given  the  name  ROSENMANN.  This  was  accomplished  by  using 
KERMIT  in  the  SE's  communications  program.  Afterwards  the  VAX 
file,  ROSENMANN,  was  edited  to  remove  blank  lines  that  were 
created  in  the  transfer.  This  data  will  be  considered  for  the  next 

version  of  the  database. 


4.7  Updates  to  the  Propram  SELECT 


There  have  been  several  changes  to  the  program  SELECT 
since  the  version  1  distribution.  In  particular  we  have  tested  in- 

house  several  different  schemes  to  make  the  program  faster.  They 
are  based  on  eliminating  as  many  of  the  code  conversions  as  possible. 
Input  to  the  program  has  been  made  easier  and  the  program  can  be 
run  in  the  interactive  mode  to  generate  a  data  file  for  the  batch 
mode.  This  work  is  aimed  at  making  the  use  of  the  database  as 
simple  as  possible. 

Other  changes  were  made  that  deal  with  problems 

associated  with  a  growing  database.  Several  times  the  array  space 

had  to  be  increased  due  to  the  addition  of  new  vibrational  bands. 

Because  the  look  up  procedure  relies  on  these  numbers,  each  time 

the  array  AVIB  is  increased  all  the  positioning  keys  need  to  be 

changed.  After  doing  this  a  few  times,  we  added  25  extra  locations 

for  the  vibrational  bands  to  each  molecule.  This  allows  us  to  add 

new  vibrational  bands  without  having  to  change  all  the  positioning 
keys. 
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5.0  PARTITION  SUMS  FOR  THE  DATABASE  SPECIES 


The  total  internal  partition  sum,  Q(T),  is  needed  to  extract 
the  transition  probability,  Rjf,  from  the  line  strength,  Sjf  (see 
equation  1).  It  is  also  used  for  converting  line  strengths  from  one 
temperature  to  another.  During  this  contract  we  started  a  program 
to  calculate  accurate  partition  sums  for  the  species  on  the  database 
and  to  consider  the  temperature  dependence  of  Q(T).  Furthermore, 
the  temperature  dependence  was  represented  by  an  analytical  form 
allowing  fast  and  accurate  calculations  of  Q(T)  without  interpolation. 
We  set  a  criterion  of  computing  Q(T)  to  better  than  \%  and  fitting 
Q(T)  vs.  T  to  better  than  1%.  For  the  molecules  considered  this  has 
been  accomplished. 


The  total  internal  partition  sum  is  given  by 


Q(T)=  Xgi-e  ‘f 


'kT 


(20) 


where  the  sum  is  over  all  states  i  of  the  molecule,  gj  is  the 

degeneracy  of  the  level  i  and  includes  only  the  nuclear  spin 
degeneracies  and  degeneracies  not  arising  from  K  or  M,  k  is  the 

Boltzmann  constant,  ej  is  the  ro-vibrational  energy  of  the  level  i,  and 

T  is  the  temperature  in  Kelvin.  We  have  evaluated  equation  20  by 

direct  summation  over  the  states  i.  We  have  studied  the  cutoff  in  i 
with  respect  to  the  accuracy  of  the  partition  sum.  What  is  needed  to 
evaluate  Q(T)  is  the  energy  levels  of  the  system,  Ej,  and  the  correct 
degeneracy  factor,  gj. 

The  degeneracy  factor  is  composed  of  several  parts. 
There  is  the  space  degeneracy  factor  (2J+1)  times  the  total  nuclear 
degeneracy  factor  of  each  state.  This  nuclear  factor  is  given  by  the 


nuclear  spin  weight  of  each  level  which  is  different  for  symmetric 
and  nonsymmetric  nuclei.  For  the  symmetric  nuclei,  the  nuclear  and 
rotational  wavefunctions  couple  to  give  a  nuclear  spin  degeneracy 
factor  (these  were  discussed  in  section  1.0).  These  are  for  example 
the  3:1  ratio  for  odd.even  states  of  H2  or  H2O.  The  non-symmetrical 
nuclei  give  rise  to  a  nuclear  spin  factor  gn  =  11(21+1),  where  I  is  the 
spin  of  each  nucleus  in  question.  In  general  one  has 
(21 1  +  l  )*(2l2+ 1  W2I3  +  I  )•••,  one  term  for  each  unsymmetrical  nuclei. 
This  factor  is  usually  omitted  because  in  most  applications  it  cancels. 
This  arises  because  for  most  applications  a  ratio  of  partition  sums  is 
used  and  the  unsymmetrical  term  is  independent  of  the  particular 
state.  In  our  early  test  calculations  we  did  not  include  the 
unsymmetrical  nuclear  spin  terms.  Once  we  were  confident  with  our 
calculations,  however,  all  terms  were  included  explicitly.  The 
expression  for  the  partition  sum  contains  the  non-symmetric  factor 
outside  of  the  summation  and  the  space  degeneracy  and  nuclear- 
rotational  wavefunction  coupling  degeneracy  inside, 


Q(T)  =  gn  X,2Ji+  11  8i 


(21) 


For  applications  where  the  factor  gn  is  needed,  it  is  important  that 
the  term  used  in  the  calculation  is  the  one  used.  Table  16  gives  the 
non-symmetrical  nuclei  contributions  for  each  species  on  the 
database  that  we  have  considered  so  far. 


With  all  the  statistics  at  hand,  only  the  energy  levels  of 
the  systems  are  needed  to  complete  the  calculations.  The  energy 
levels  that  are  needed  can  be  obtained  from  the  database.  Care  was 
exercised  to  extract  all  levels  but  only  for  unique  quantum  states  of 
the  molecules.  To  obtain  accurate  sums  no  duplications  of  levels  can 
occur.  The  approach  was  to  study  the  band  sums  for  the  possible 
lower  vibrational  states  of  each  molecule  to  be  considered.  For  each 
of  the  vibrational  states,  the  energy  level  corresponding  to  each 
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Table  16.  Nuclear  Spin  Degeneracy  Factors  for 
Non-Symmetrical  Nuclei. 


AFGL 

Isotope  gj  =  n ( 2 1  j  +  1 ) 

Molecule  Code 


H20  161  1 

181  1 

171  6 

162  6 

CCh  626  1 

636  2 

628  1 

627  6 

638  2 

637  12 

828  1 

728  6 

03  666  1 

668  1 

686  1 

N20  446  9 

456  6 

546  6 

448  9 

447  54 

CO  26  I 

36  2 

28  1 

27  6 

CH4  211  1 

311  2 

212  3 

N02  646  3 


unique  rotational  state  is  taken  and  written  to  a  file  for  that 
molecule.  Because  of  the  structure  of  the  database,  this  created  two 
groups,  one  for  the  molecules  that  have  the  rotational  quantum 
numbers  in  QLOCLO,  e.g.  H2O,  HNO3,  etc.,  and  the  other  for  molecules 
that  have  branch  information  and  quantum  numbers  in  QLOCLO,  e.g. 
CO2,  OH,  etc.  The  problem  is  further  complicated  by  the  presence  of 
hyperfine  structure  for  some  of  the  molecules.  Thus  rotational 
groups  1,  3,  and  4  on  the  database  (Table  1)  store  quantum  numbers 
so  QLOCLO  values  can  be  compared  and  stored  since  they  uniquely 
define  the  state.  For  group  2  molecules  (branch  information  and 
integer  J")  J"  and  F”  were  extracted  and  written  to  a  character  string 
and  these  variables  were  compared  and  saved  since  they  are  unique. 
For  group  5  (oxygen  molecule)  we  stored  J",  N",  and  F"  in  the 

character  string  and  for  group  6  RJ"  and  F"  were  stored  and  similar 
comparisons  done.  A  single  program  was  written  to  perform  this 
task  and  run  on  the  database.  This  generated  files  of  rotational  state 
information  and  lower  state  energies  for  each  vibrational  state  of 
each  isotopic  species  of  each  molecule.  The  files  were  printed  and  a 
copy  delivered  to  L.  Rothman.  These  files  revealed  some  format 
inconsistencies  that  will  be  addressed  on  the  next  version  of  the 
database.  This  resulted  in  energy  levels  to  use  to  compute  the 

partition  sums.  Preliminary  calculations  were  done  for  H2O  and  CO2 
at  300K  and  compared  with  published  values.  At  300K  for  12(M6o2 
we  found  Q  =  291.1136  compared  with  the  value  of  Gray  and 
Young29  of  291.0471.  For  H2^0,  we  found  Q(300)  =  154.0377 
compared  with  the  JPL  catalogue^  value  of  178.115.  After 
discussions  with  H.  Pickett  of  JPL,  we  agreed  our  value  was  the  more 
correct  one. 

We  next  considered  what  cutoff  in  energy  was  required 

to  produce  the  partition  sum  accurate  to  better  than  1%.  This  was 

studied  by  plotting  the  computed  Q(T)  vs.  T  for  H2O  and  CO2  at 
several  temperatures  and  seeing  is  the  partition  sum  has  converged. 
Figures  5-16  are  plots  for  the  main  isotopic  species  of  H2O  (Figures 
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5-10)  and  for  CO2  (Figures  11-13).  The  figures  show  that  the 
partition  sums  have  converged.  When  similar  plots  (Figures  14-16) 
for  the  lesser  isotopic  species  of  CO2  were  done  it  became  clear  that 
the  calculations  did  not  converge. 

This  suggested  that  we  did  not  have  the  "complete"  set  of 
states  available  to  us  from  the  energy  levels  extracted  from  the 
database.  This  was  apparent  for  the  lesser  isotopic  species  but  on 
closer  inspection  it  was  also  true  for  the  principal  species.  Our 
analysis  indicated  the  partition  sum  could  be  very  accurately 
determined  by  including  all  levels  up  to  5000  wavenumbers.  This 
ensures  accurate  values  for  Q(T)  up  to  temperatures  of  500K.  To 
show  this  we  extracted  the  energy  levels  as  a  function  of  rotational 
quantum  index  (J(J+1)  +Ka-Kc+1)  for  the  principal  species  of  water 
vapor  and  ozone.  These  are  displayed  in  Tables  17  and  18.  The 
zeros  in  the  tables  are  when  the  database  has  no  energy  value  for 
the  state.  This  clearly  showed  that  we  could  not  use  the  energy 
levels  from  the  database  only  and  expect  partition  sum  to  better 
than  [%. 

We  discussed  different  schemes  to  produce  the  needed 
energy  levels.  The  approximate  schemes  that  we  tried  could  not 
produce  the  energy  levels  with  the  accuracy  that  was  needed.  At 
this  point  we  abandoned  the  use  of  approximate  schemes  and 
decided  to  calculate  all  energy  levels  from  molecular  Hamiltonians. 

The  other  aspect  of  these  calculations  is  how  to  represent 
the  partition  sums  as  a  function  of  temperature.  The  method 
commonly  used  is  to  calculate  Q(T)  for  a  set  number  of  temperatures 
and  to  linearly  interpolate  on  a  In { Q(T) }  vs.  In { T }  basis.  This  is  not  an 
elegant  method  to  get  Q(T)  at  arbitrary  temperatures  but  is  certainly 
applicable  if  all  else  fails.  We  have  sought  simple  analytical 
expressions  that  can  represent  Q(T)  in  terms  of  a  small  number  of 
parameters.  Our  first  tests  were  on  a  simple  two  parameter  power 
law,  Q(T)  =  aT^.  This  worked  very  well  for  the  principal  species  of 
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water  vapor  from  200-350K,  0.04%  maximum  error  in  reproducing 
Q(T).  This  form  failed  when  applied  to  other  isotopic  species  and  to 
other  molecules. 

A  more  complex  analytical  form  based  on  the  work  of 
McDowell31  was  tried  next.  From  McDowell,  the  analytical  form  of 
the  rotational  partition  sum  for  a  spherical  top  molecule  is  given  by 


Qr  =  aY/2 


P/ 


-3 


1  + 


15D/ 


'43B/ 


(22) 


where  p  =  hcB/kT  (B  =  rotational  constant)  and  a*  =  (2Iy  +  l)n/a  (Iy  = 
nuclear  spin  of  the  Y  nuclei,  o  =  classical  symmetry  number)  and  D  = 
centrifugal  distortion  constant.  In  this  expression  a,  p,  D  and  B  were 
treated  as  adjustable  constants.  This  proved  to  give  good  results. 

We  then  tried  a  simple  polynomial  expression  in 
temperature  with  four  adjustable  constants, 


2  3 

Q(T)  =  a+bT+cT  +dT  ,  (23) 

This  form  for  the  partition  sum  proved  to  be  the  simplest  for 
recalculation  of  Q(T)  and  it  gave  the  most  accurate  results.  This  was 
the  form  that  we  chose  to  represent  the  partition  sums. 

We  then  fixed  the  temperature  range  of  study  to  consider 
from  70  to  440K  in  5K  steps.  The  energy  levels  for  each  molecule 
were  calculated  via  the  appropriate  Hamiltonian  for  each  vibrational 
state  and  the  partition  sums  done.  The  partition  sums  were  then 
f i tted ^ 2  by  the  polynomial  expression  by  first  doing  a  least  squares 
fit.  The  results  of  the  least  squares  fit  was  used  as  input  to  a 
minimax  (Simplex)  procedure  to  the  eight  power  (the  power  of  eight 
was  chosen  by  increasing  the  exponent  until  no  further  change  in  the 
results  was  observed).  The  standard  deviation  of  the  least  squares 
fit  was  used  as  a  step  size  in  the  Simplex  minimization.  This  yields  a 
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set  of  coefficients;  a,  b,  c,  d,  for  each  isotopic  species  of  each  molecule 
studied.  Below  we  discuss  each  molecule  separately. 


5.1  Water  Vapor 

For  water  vapor,  there  are  some  levels  missing  for  the 
principal  isotopic  species  (161  AFGL  code)  and  many  levels  missing 
for  the  other  species  181,  171  and  162.  By  using  the  index,  J  +  Ka  - 
Kc  +  1,  we  have  filled  an  array  of  energy  levels  for  each  isotopic 
species.  This  allows  us  to  determine  which  levels  are  missing  for 
each  isotopic  species.  We  then  took  the  Watson  Hamiltonian 
constants  from  C.  Camy-Peyret  et.  al. 33-34  for  the  161,  181  and  171 
isotopic  species  and  the  Hamiltonian  constants  of  Toth  et.  al.35  for  the 
HDO  species  (162)  (the  constants  are  given  in  Table  19)  and  used 
these  to  generate  all  the  missing  levels.  In  order  to  do  this  the 
asymmetric  rotor  energy  program  had  to  be  made  operable.  After 
some  work  this  was  done  and  test  calculations  were  made  to  check 
the  integrity  of  the  program.  Calculations  were  made  up  to  a  cutoff 
in  energy  which  relates  to  the  maximum  J  level  for  a  vibrational 
state.  This  gave  us  a  near  complete  set  of  energy  levels  up  to  5000 
cm'l.  Because  the  energy  levels  from  the  database  come  from  the 
work  of  Camy-Peyret  et  al.,  the  final  energies  are  a  consistent  set. 

The  partition  sums  were  then  calculated  using  the 
"complete"  set  of  energy  levels.  The  calculations  were  done  in  part 
on  the  AFGL  VAX  and  on  an  ULCAR  IBM-AT-PC.  In  Figure  17  the 
convergence  of  the  partition  sum  at  several  temperatures  is  shown. 
The  partition  sum  has  leveled  off  at  around  2500  cm-1.  Error 

analysis  was  also  done  by  shifting  the  energy  levels  by  the  maximum 
error  reported  in  the  literature  (0.2%),  recalculating  the  partition 
sums  and  then  forming  the  difference  of  the  two  values.  The  results 


Flaud  arid  Camy-Peyret  ground  state  161  Watson  constants 


0 . 00G00000 

E+00 

Vibrational  energy 

27 . 880678 

E+00 

A 

14. 521689 

E+00 

B 

9.277459 

E+00 

C 

3 .25199 

E  -  02 

DLK 

-5.7655 

E-  03 

DLJK 

1 . 24894 

E-03 

DLJ 

1 . 3007 

E  -  03 

dsK 

5.0838 

E  -  04 

dsJ 

1.2261 

E  -  04 

HLK 

-1.3700 

E  -  05 

HLKJ 

-1.921 

E  -  06 

HUK 

4 . 5061 

E-07 

HU 

2 . 3958 

E  -  0  5 

hsK 

-4  .21 

E-07 

hsJK 

2.2568 

E-07 

hsJ 

-6.496 

E-07 

ALK 

1 .  29150 

E-07 

ALKKJ 

-3.681 

E  -  08 

AUK 

0.0 

E  -  00 

AU 

8 . 66000 

F.-10 

AUJK 

-  1  .  5414 

E-07 

AlsK 

6  .  1  7 

E-  10 

AlsJK 

0 . 0 

E  -  00 

AlsKJ 

0.0 

E-  00 

AlsJ 

2  .  399 

E  -  0  9 

PK 

6  .  8 1 

E-l? 

PKKJ 

1  .  704 

E  - 10 

psK 

-6.49 

E- 12 

QK 

7.65 

E  - 1 5 

RK 

■;  :+.  s'  n  Hamiltonian  Constants  for 
Me  o t > r : c i o s  of  I  i  „ ' 


ground  state  181  Watson  constants 

Vibrational  energy 
A 
B 
C 

DLK 

DLJK 

DLJ 

dsK 

dsJ 

HLK 

HLKJ 

HLJK 

HLJ 

hsK 

hsJK 

hsJ 

ALK 

ALKKJ 

ALJK 

AU 

A1JJK 

AlsK 

Als.JK 

Al  ■  KJ 

Al  sJ 

PK 

PKKJ 

psK 

QK 

RK 


Flaud  and  Camy-Peyret  ground  state  171  Watson  constants 


0.00000000 

E+00 

27.695220 

E+00 

14 . 521875 

E+00 

9 .256840 

E+00 

3.20133 

£-02 

-5.71998 

E-03 

1.25345 

E-03 

1.2885 

E-03 

5.0841 

E  -  04 

1 . 2040 

E-  04 

-1.6324 

E-  05 

-1.378 

E-06 

5.3645 

E-  07 

3.1329 

E  -  05 

-7.396 

E-07 

2.7281 

E-  07 

-6.882 

E-07 

2.4318 

E-07 

-6.947 

E-08 

-2.759 

E  - 10 

0.00000 

E  00 

-2. 8704 

E-07 

0.00 

E  00 

-7.76 

E-09 

-1.533 

E  - 10 

2.031 

E-09 

0.00 

E  00 

4.613 

E  - 10 

-3.964 

E- 12 

0.00 

E- 15 

Vibrational  energy 
A 
B 
C 

DLK 

DUK 

DU 

dsK 

dsJ 

HLK 

HLKJ 

HUK 

HU 

hsK 

hsJK 

hsJ 

ALK 

ALKKJ 

AUK 

AU 

AUJK 

AlsK 

AlsJK 

AlsKJ 

AlsJ 

PK 

PKKJ 

psK 

QK 

RK 


Table  11.  (Continued) 
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Toth's  Watson  constants  for  HDO 


0.00000000 

E+00 

Vibrational  energy 

23.413842 

E+00 

A 

9 .103323 

E+00 

B 

6.406295 

E+00 

C 

1.25376 

E-02 

DLK 

1 . 1432 

E-03 

DLJK 

3.60562 

E-  04 

DLJ 

2.0814 

E-03 

dsK 

1.2143 

E  -  04 

dsJ 

4 .5589 

E-05 

HLK 

-8.702 

E-06 

HLKJ 

2.401 

E  -  06 

HLJK 

3.513 

E-08 

HU 

1.482 

E-05 

hsK 

1.001 

E-06 

hsJK 

2.112 

E-08 

hsJ 

-1.055 

E-07 

ALK 

0.00000 

E  00 

ALKKJ 

0.00000 

E  00 

AUK 

0.000 

E-08 

ALK 

0.00000 

E  00 

AU.1K 

0.00000 

E  00 

AlsK 

0.00000 

E  00 

AlsJK 

0  00000 

E  00 

AlsKJ 

1.09000 

E- 11 

AlsJ 

0.00000 

E  00 

PK 

0.00000 

E  00 

PKKJ 

0.00000 

E  00 

psK 

0.00000 

E  00 

QK 

0.00000 

E  00 

RK 

Tab  1  e 
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7 
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! 
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6 

of  the  error  analysis  show  the  partition  sums  accurate  to  -0.3%, 
which  is  certainly  in  the  error  range  we  desired. 

The  polynomial  fits  were  done  for  the  four  isotopic 
variants  of  water  vapor.  Figures  18  and  19  display  the  %  error  |(Q(T) 
-  polynomial)*  100/Q(T)]  of  the  fits  for  and  for  HDO.  For  the 

principal  species  the  error  remains  below  0.25%,  for  HDO  it  remains 
below  0.02%,  confirming  the  validity  of  the  polynomial  fit.  Thus  we 
are  confident  that  the  partition  sums  are  accurate  to  better  than  1%. 


5.2  Carbon  Dioxide 


For  carbon  dioxide,  all  energy  levels  upto  5000  cm-' 
were  calculated  via  L.  Rothman’s  C02-Generator  program.  This 
program  used  the  most  recent  set  of  molecular  constants  for  the 
isotopic  species  13(M6o2,  and  12(^180.  The  error 

analysis  used  a  shift  in  the  energy  levels  of  0.2%,  this  is  greater  than 
the  actual  error  for  most  of  the  lines  but  provides  an  upper  limit. 
With  this  shift  in  energy  the  partition  sums  did  not  change  by  more 
than  0.3%.  The  partition  sums  were  evaluated  at  the  specified 
temperatures  and  then  the  results  fitted  by  the  polynomial.  The 
resulting  error  of  the  fitting  is  demonstrated  in  Figure  20  for  the 
principal  isotopic  species  of  C02-  Note  the  fitting  error  remains  less 
than  -0.02%  again  confirming  the  polynomial  fit  and  total  internal 
partition  sums  accurate  to  better  than  1%. 


5.3  Ozone 


Ozone  had  many  energy  levels  missing  as  well  as  having 
many  states  coming  from  older  less  reliable  sources.  Because  better 
values  are  available,  especially  for  the  isotopic  species,  we  decided  to 


DI*D/(lBIUIOUAl°aOIBD)*OOI 

10XX&  % 


,  Error  | 0 (T ) -Polynomial /QT 


vs . 


OIED/(lT3T\«O\xApd-OpD)*00I 

*°«3% 


(l)0/(ieiwouAiod-(i)0)*‘OOl 
0Duajaiiia  % 


Error  | Q (T ) -Polynomial /QT |  vs.  T  for  CO 


do  the  ozone  partition  sums  from  molecular  constants  entirely,  i.e. 
the  levels  from  HITRAN  will  not  be  used  at  all.  The  references  for 
the  ozone  Watson  Hamiltonian  constants  are  Flaud  et.  al. 36-41  These 
constants  are  (see  Table  20)  for  the  principal  species  (666)  for  the 
states  000,  001,  100,  110,  Oil,  010,  020,  002,  101  and  200.  For  me 
686  and  668  species  the  constants  include  the  000,  010,  001  and  100 
vibrational  states.  Note  for  the  principal  isotopic  species  (666)  there 
are  ten  vibrational  levels  and  we  tested  the  often  used 
approximation 


Q(T)  =  QV(T)  •  Qr(T)  .  (24, 

used  for  these  levels.  Many  thousands  of  levels  were  calculated  by 
the  ASMROT  program  and  formatted  for  the  partition  sum  program. 
Again  the  5000  cm'1  cutoff  was  enforced  to  produce  the  final  set  of 
"complete"  energy  levels.  The  resulting  partition  sums  and  the  error 
analysis  are  given  in  Tables  21  and  22.  The  error  analysis  shows  the 
partition  sums  being  good  to  better  than  0.4%  error. 

The  test  of  the  approximation  of  Qvib  times  Qrot  being 
equal  to  Qtot  appears  very  good  for  ozone.  The  difference  between 
the  actual  partition  sum  and  the  approximation  was  never  greater 
than  0.16%  error  from  100  to  400  K  (see  Table  23). 


5.4  N?Q 


For  N2O  we  sorted  the  energy  levels  from  the  database 
and  found  many  levels  missing  for  J  >  70.  These  energies  are  in  the 
range  2500  cm'1  and  must  be  accounted  for  if  the  calculations  are  to 
be  of  the  required  accuracy.  The  vibrational  bands  of  importance 
are: 


o2 


Flaud  et.  al .  ozone  (666)  ground  state  Watson  constants 


0.00000000 

E+00 

3.553666502 

E+00 

0 . 445283207 

E+00 

0.394751807 

E+00 

0.2116872 

E-  03 

0 . 184846 

E-05 

0.454169 

E-06 

0.323240 

E-05 

0.697903 

E-07 

0.39896 

E  -  07 

0.18592 

E  -  08 

0.646 

E  - 11 

0.321 

E- 12 

0.2293 

E-08 

0.757 

E- 11 

0.1755 

E- 12 

0.1177 

E- 10 

0.273 

E- 12 

0.0 

E-00 

0.299 

E- 17 

0.325 

E- 15 

0.181 

E- 12 

0.0 

E-00 

0.0 

E  00 

0.0 

E  00 

C.223 

E  - 14 

0.189 

E- 15 

0.0 

E-00 

0.0 

E-00 

0.0 

E-00 

Vibrational  energy 
A 
B 
C 

DLK 

DUK 

DU 

dsK 

dsJ 

HLK 

HLKJ 

HUK 

HU 

hsK 

hsJK 

hsJ 

ALK 

ALKKJ 

AUK 

AU 

AUJK 

AlsK 

AlsJK 

AlsKJ 

AlsJ 

PK 

PKKKJ 

psK 

QK 

RK 


20. 


Watson  Hamiltonian  Constants  for  C, 
from  Fl-vrl  et  al. 


Flaud  et.  al .  ozone  (666)  010  state  Watson  constants 


700.931056 

E+00 

3.607092966 

E+00 

0.444021809 

E  +  00 

0.392439113 

E+00 

0.2327460 

E  -  03 

-0.1790036 

E-05 

0.457106 

E  -  06 

0.388935 

E-05 

0.692601 

E  -  07 

0.48866 

E  -  07 

-C. 24209 

E  -  08 

0.309 

E- 11 

0.31367 

E- 12 

0.31968 

E-08 

-0.1155 

E- 10 

0.1868 

E- 12 

-0.1582 

E- 10 

0.4795 

E  - 12 

0.0 

E-00 

0.0 

E- 17 

0.0 

E- 15 

0 .656 

E- 12 

0.0 

E-00 

0.383 

E  - 13 

A  A 

U  .  \J 

E  00 

0.498 

E  - 14 

0.0 

E- 15 

0.0 

E-00 

0.0 

E-00 

0.0 

E-00 

Vibrational  energy 
A 
B 
C 

DLK 

DDK 

DLJ 

dsK 

dsJ 

HLK 

HLKJ 

KLJK 

HU 

hsK 

hsJK 

hsJ 

ALK 

ALKKJ 

AUK 

AL) 

AUJK 

AlsK 

AlsJK 

AlsKJ 

AlsJ 

PK 

PKKKJ 

?sK 

QK 

RK 


Tabic  20.  (  Tent  i:v;c  ) 
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Flaud  ec.  al .  ozone  (666)  001  state  Watson  constants 


1042.08398 

E+00 

3 . 50055300 

E+00 

0.44129635 

E+00 

0.390998275 

E+00 

0 . 208072 

E  -  03 

0.16977 

E  -  05 

0.456697 

E-06 

0.29609 

E-  05 

0.759513 

E-07 

0.39036 

E-07 

0.1802 

E-08 

0.160 

E- 11 

0.1704 

E- 12 

0.248 

E-08 

0.1142 

E- 10 

0.2355 

E- 12 

0.1141 

E  - 10 

0.317 

E- 12 

0.0 

E-  00 

0. 2994 

E- 17 

0.3251 

E- 15 

0.21 

E  - 12 

C.O 

E-00 

0.0 

E  - 13 

0.0 

E  00 

0.333 

E  - 14 

0.1895 

E- 15 

0.0 

E-00 

0.0 

E-00 

0.0 

E-00 

Vibrational  energy 
A 
B 
C 

DLK 

DUK 

DU 

dsK 

dsJ 

HLK 

HLKJ 

HUK 

HU 

hsK 

hsJK 

hs  J 

ALK 

ALKKJ 

AUK 

AU 

AUJK 

AlsK 

AlsJK 

AlsKJ 

AlsJ 

PK 

PKKKJ 

psK 

QK 

RK 


l  a: 


(Continue 'I ) 


Flaud  et.  al .  ozone  (666)  100  state  Watson  constants 


1103 . 13729 

E+00 

Vibrational  energy 

3.55669374 

E+00 

A 

0.442736468 

E+00 

E 

0 . 392568222 

E+00 

C 

0.217495 

E-03 

DLK 

-0.21569 

E-  05 

DUK 

0.461481 

E  -  06 

-  DLJ 

0.35570 

E-05 

dsK 

0.658775 

E-07 

risJ 

0.41897 

E  -  07 

hu: 

-0.2041 

E  -  08 

HLKJ 

-0.312 

E  - 1 0 

HUK 

0 . 424 

E  - 12 

HU 

0.98 

E  -  09 

hsK 

-0.80 

E- 11 

hsJK 

0.916 

E- 13 

hsJ 

-0.12189 

E  - 10 

AEK 

0.309 

E  - 12 

ALKKJ 

0.0 

E-00 

AUK 

0.2994 

E- 17 

au 

-0.3251 

E  - 15 

A  v1  K 

0.15 

E- 12 

Al  sK 

o 

o 

E-00 

AlsJK 

o 

o 

E- 13 

AlsKJ 

o 

o 

E  00 

Al  s  J 

0.333 

E  - 14 

PK 

0.1895 

E- 15 

PKKKJ 

o 

o 

E-00 

psK 

0.0 

E-00 

QK 

0.0 

E-00 

RK 

(:  -  i  r.u-  : 


Flaud  et.  al .  ozone  '666)  020  state  Vatson  constants 


1399.272626 

E+00 

3.6624062 

E+00 

0.44274468 

E+00 

0.39007214 

E+00 

0.256077 

E-03 

0.17212 

E-05 

0.46006 

E-06 

0.45898 

E-05 

0.68637 

E-  07 

0 .610 

E-07 

0.2542 

E-08 

0.186 

E  - 10 

0.118 

E- 12 

0.478 

E-08 

0.268 

E  - 10 

0.238 

E- 12 

0.333 

E- 10 

0.953 

E- 12 

0.0 

E-00 

0.0 

E- 17 

O.C 

E- 15 

0.0 

E- 12 

0.0 

E-00 

0.0 

E- 13 

0.0 

E  00 

C  .  40 

E-13 

0.0 

E-00 

0.0 

E-00 

0.0 

E-00 

0.0 

E-00 

Vibrational  energy 
A 
B 
C 

DLK 

DUK 

DLJ 

dsK 

dsJ 

HLK 

HLKJ 

HLJK 

HU 

hsK 

hsJK 

hsJ 

ALK 

ALKKJ 

AUK 

AU 

AUJK 

AlsK 

AlsJK 

AlsKJ 

AlsJ 

PK 

PKKKJ 

psK 

QK 

RK 


Table  20.  (Continued) 
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Flaud  et.  al.  ozone  (666)  002  state  Watson  constants 


2063.17411 

E+00 

3.44883 

E+00 

0.4374296 

E+00 

0.367865 

E+00 

0.20455 

E-03 

-0.1426 

E-05 

0.46727 

E-06 

0.2075 

E-05 

0.7886 

E-  07 

0.382 

E  -  07 

-0 . 1474 

E-08 

0.0 

E- 11 

0.3215 

E- 12 

0.327 

E-08 

-0.7576 

E- 11 

0.1755 

E- 12 

-0.104 

E- 10 

0.2739 

E- 12 

0.0 

E-00 

0.2994 

E- 17 

-0.2251 

E- 15 

0.1812 

E- 12 

0.0 

E-00 

0.0 

E-13 

0.0 

E  00 

0.333 

E-14 

0.1895 

E- 15 

0.0 

E-00 

0.0 

E-00 

0.0 

E-00 

Vibrational  energy 
A 
B 
C 

DLK 

DLJK 

DLJ 

dsK 

dsJ 

HLK 

HLKJ 

HUK 

HU 

hsK 

hsJK 

hsJ 

ALK 

ALKKJ 

AUK 

AU 

r.UJK 

AlsK 

AlsJK 

AlsKJ 

AlsJ 

PK 

PKKKJ 

psK 

QK 

RK 


Table  20. 


( Continued ) 


8  8 


Flaud  et.  al .  ozone  (666)  101  state  Watson  constants 


2110.78428 

E+00 

3.5019861 

E+00 

0.4385779 

E+00 

0.388266 

E+00 

0.213845 

E-  03 

0.1870 

E-05 

0.46090 

E-06 

0.4330 

E-05 

0.7482 

E  -  07 

0.4085 

E-07 

0.2173 

E-08 

0.22 

E  - 10 

0.3215 

E- 12 

0.207 

E-08 

0.7576 

E- 11 

0.1755 

E- 12 

0.1106 

E- 10 

0.2739 

E- 12 

0.0 

E-00 

0.2994 

E- 17 

0.3251 

E  - 1 5 

0.1612 

E- 12 

0.0 

E-00 

0.0 

E- 13 

C.O 

E  00 

0.333 

E-14 

0.1895 

E- 15 

C.O 

E-00 

0.0 

E-00 

0.0 

E-00 

Vibrational  energy 

A 

B 

C 

DLK 

DUK 

DU 

dsK 

dsJ 

HLK 

HLKJ 

HUK 

HU 

hsK 

hsJK 

hsJ 

ALK 

ALKKJ 

AUK 

AU 

AUJK 

AlsK 

AlsJK 

AlsKJ 

AlsJ 

PK 

PKKKJ 

psK 

QK 

RK 


Table  2  0.  (Continued) 
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Flaud  et.  al .  ozone  (666)  200  state  Watson  constants 


2195.87187 

E+00 

3.55982 

E+00 

0 . 4400867 

E+00 

0.390089 

E+00 

0.22394 

E-03 

-C . 2641 

E-  05 

0.46579 

E-06 

0.3510 

E-05 

0.6250 

E  -  0  7 

0.468 

E  -  07 

-0.2405 

E-  08 

-0.57 

E- 10 

0.3215 

E- 12 

C.  143 

E-08 

-0.7576 

E- 11 

0.1755 

E  - 12 

-0.191 

E- 10 

0.2739 

E- 12 

0.0 

E-00 

0.2994 

E-17 

-  0.3251 

E- 15 

0.1812 

E- 12 

0.0 

E-00 

0.0 

E- 13 

0.0 

E  00 

0.333 

E*  14 

0.1895 

E  - 15 

0.0 

E-00 

0.0 

E-00 

0.0 

E-00 

Vibrational  energy 

A 

B 

C 

DLK 

DUK 

DLJ 

dsK 

dsJ 

HLK 

HLKJ 

HUK 

HU 

hsK 

hsJK 

hsJ 

ALK 

ALKKJ 

AUK 

AU 

AUJK 

AlsK 

AlsJK 

AlsKJ 

AlsJ 

FK 

pkkkj 

psK 

QK 

RK 


iablc  20.  (Con-rinuod) 


Flaud  et.  al.  ozone  (666)  110  state  Vatson  constants 


1796.26187 

E+00 

3.6108121 

E+00 

0.44144616 

E+00 

0 . 39029670 

E+00 

0.239727 

E-03 

-0.20753 

E-05 

0.46689 

E-  06 

0.4321 

E-05 

0.6473 

E-07 

0.4997 

E-07 

-0.2104 

E  -  08 

-0.84 

E  - 10 

0.169 

E- 12 

-0.28 

E-08 

-0.33 

E- 10 

0.92 

E- 13 

-0.126 

E- 10 

0.2739 

E- 12 

0.0 

E-  00 

0.2994 

E- 17 

-0.3251 

E- 15 

0.1612 

E- 12 

0.0 

E-00 

0.0 

E- 13 

0.0 

E  00 

0.333 

E- 14 

0.1895 

E- 15 

0.0 

E-00 

0.0 

E-00 

0.0 

E-00 

Vibrational  energy 
A 
B 
C 

DLK 

DUK 

DU 

dsK 

dsJ 

HLK 

HLKJ 

HUK 

HU 

hsK 

hsJK 

hsJ 

ALK 

ALKKJ 

AUK 

AU 

AUJK 

AlsK 

AlsJK 

AlsKJ 

AlsJ 

PK 

PKKKJ 

psK 

QK 

RK 


Table  21.  (Continued) 
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Flaud  et.  al .  ozone  (666)  011  state  Watson  constants 


1726.52249 

E+00 

Vibrational  energy 

3.5523057 

E+00 

A 

0.43989712 

E+00 

B 

0.38850339 

E  +  00 

C 

0.22865 

E-03 

DLK 

-0.1596 

E-05 

DUK 

0.45985 

E-06 

DU 

0.3495 

E-05 

asK 

0.7610 

E-07 

dsJ 

0.464 

E-  07 

HLK 

-0.1876 

E-08 

HLKJ  | 

-0.311 

E- 10 

HUK 

0.114 

E- 12 

HU 

0.0 

E-08 

hsK 

0.0 

E  - 11 

hsJK 

0.200 

E- 12 

hsJ 

-0.108 

E- 10 

ALK 

0.2739 

E- 12 

ALKKJ 

0.0 

E-00 

AUK 

0.2994 

E- 17 

AU 

-  0 . 32  51 

E- 15 

AUJK 

0.1812 

E- 12 

AlsK 

0.0 

E-00 

AlsJK 

0.0 

E- 13 

AlsKJ 

0.0 

E  00 

AlsJ 

0.333 

E- 14 

PK 

0.1895 

E- 15 

PKKKJ 

0.0 

E-00 

psK 

0.0 

E-00 

QK 

0.0 

E-00 

RK 

31 


Flaud  et.  al .  ozone  (686)  ground  state  Watson  constants 


0.00000000 

E+00 

Vibrational  energy 

3.29049900 

E+00 

A 

0.445399229 

E+00 

B 

0.391329662 

E+00 

C 

0.1808811 

E-03 

DLK 

0.13179 

E  -  05 

DUK 

0.44784 

E-06 

DU 

0.31155 

E-05 

dsK 

0.72905 

E-07 

dsJ 

0.29550 

E  -  07 

HLK 

0.1475 

E-08 

HLKJ 

0.121 

E- 10 

HUK 

0.439 

E- 12 

HU 

0.161 

E-08 

hsK 

0.11 

E- 11 

hsJK 

0.160 

E- 12 

hsJ 

o 

o 

E-00 

ALK 

o 

o 

E-00 

ALKKJ 

0.0 

E-00 

AUK 

0.0 

E-00 

ALK 

o 

o 

E-00 

AUJK 

0.0 

E-00 

AlsK 

0.0 

E-00 

AlsJK 

0.0 

E  00 

AlsKJ 

0.0 

E  00 

AlsJ 

0.0 

E-00 

PK 

0.0 

E-00 

PKKJ 

0.0 

E-00 

psK 

0.0 

E-00 

QK 

0.0 

E-00 

RK 

Table  2C.  (Continued) 
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Flaud  et.  al .  ozone  (686)  010  state  Watson  constants 


693.30570 

E+00 

3 . 34155121 

E+00 

0.444188088 

E+00 

0.388900751 

E+00 

0.1990522 

E-03 

0.12424 

E  -  05 

0.450826 

E  -  06 

0.37196 

E-05 

0.72243 

E-07 

0.37749 

E  -  07 

0.1709 

E  -  08 

0.25 

E  - 11 

0.332 

E- 12 

0.262 

E-08 

0.55 

E- 11 

0.173 

E- 12 

0.674 

E- 11 

0.0 

E-00 

0.0 

E-00 

0.0 

E-00 

c.  o 

E-00 

0.0 

E-00 

0.0 

E-00 

0.0 

E  00 

0.0 

E  00 

0.0 

E-00 

0.0 

E-00 

0.0 

E-00 

0.0 

E-00 

0.0 

E-00 

Vibrational  energy 
A 
B 
C 

DLK 

D’JK 

DIJ 

dsK 

dsJ 

hu: 

HLKJ 

HUK 

HU 

hsK 

hsJK 

hsJ 

ALK 

ALKKJ 

AUK 

ALK 

AUJK 

AlsK 

Al  s  JK 

AlsKJ 

AlsJ 

PK 

PKKJ 

psK 

QK 

RK 


Table-  2  0 


(Continued) 
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Flaud  et.  al.  ozone  (686)  100  state  Watson  constants 


1074.30756 

E+00 

3.290456 

E+00 

0.4428737 

E+00 

0.386777 

E+00 

0.181 

E-03 

-0.132 

E-05 

0.448 

E-06 

0.311 

E-05 

0.729 

E-07 

0.313 

E-07 

-0.158 

E-08 

-0.300 

E  - 11 

0.327 

E- 12 

0.217 

E-08 

-0.967 

E- 11 

0.167 

E- 12 

0.0 

E-00 

0.0 

E-00 

0.0 

E-00 

0.0 

E-00 

0.0 

E-00 

0.0 

E-00 

0.0 

E-00 

0.0 

E  00 

0.0 

E  00 

O.C 

E-00 

0.0 

E-00 

0.0 

E-00 

0.0 

E-00 

0.0 

E-00 

Vibrational  energy 
A 
B 
C 

DLK 

DUK 

DLJ 

dsK 

dsJ 

HLK 

HLKJ 

HUK 

HLJ 

hsK 

hsJK 

hsJ 

ALK 

ALKKJ 

AUK 

ALK 

AUJK 

AlsK 

AlsJK 

AlsKJ 

AlsJ 

PK 

PKKJ 

psK 

QK 

RK 


Table  20.  (Continued) 
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Flaud  et.  ai .  ozone  (686)  001  state  Vatson  constants 


1008 .45277 

E+00 

Vibrational  energy 

3.243579 

E+00 

A 

0.4415969 

E+00 

E 

0.386332 

E+00 

r 

0.17836 

E-03 

DLK 

-0.1387 

E  -  05 

DUK 

0.45380 

E  -  06 

DU 

0.2604 

E-05 

dsK 

0.7713 

E-  07 

dsJ 

0  3062 

E-07 

HLK 

-0.99 

E-  09 

HLKJ 

-0.300 

E  - 11 

HUK 

0.327 

E- 12 

HU 

0.217 

E-08 

hsK 

-0.967 

E  - 11 

hsJK 

0.187 

E- 12 

hsJ 

0.0 

E-00 

ALK 

0.0 

E-  00 

ALKKJ 

0.0 

E-00 

AUK 

0.0 

E-00 

ALK 

0.0 

E-00 

AUJK 

o 

O 

E-00 

AlsK 

0.0 

E-00 

AlsJK 

0.0 

E  00 

AlsKJ 

0.0 

E  00 

AlsJ 

0.0 

E-00 

PK 

0.0 

E-00 

PKKJ 

0.0 

E-00 

psK 

0.0 

E-00 

QK 

0.0 

E-00 

RK 

Tatle  20.  (  Cor/. ) 


06 


Flaud  et.  al.  ozone  (668)  ground  state  Vatson  constants 


0.00000000 

E+00 

Vibrational  energy 

3.48818507 

E+00 

A 

0.420008328 

E+00 

B 

0 . 374008935 

E+00 

C 

0.2038645 

E  -03 

DLK 

-0.190229 

E-05 

DUK 

0.406917 

E-06 

DU 

0.29205 

E-05 

dsK 

0.606059 

E-  07 

dsJ 

0.35200 

E-07 

HLK 

-0.1589 

E-  08 

HLKJ 

-0.140 

E  - 10 

HUK 

0.300 

E- 12 

HU 

0.159 

E  -  08 

hsK 

0.20 

E  - 1 1 

hsJK 

0.112 

E- 12 

hsJ 

c.o 

E-  00 

ALK 

c.o 

E  -  00 

ALKKJ 

0.0 

E-00 

AUK 

0.0 

E-00 

ALK 

0  0 

E-00 

AUJK 

o 

o 

E-00 

AlsK 

o 

o 

E-00 

AlsJK 

0.0 

E  00 

AlsKJ 

0.0 

E  00 

Al  s  J 

0.0 

E-00 

PK 

0.0 

E-00 

PKKJ 

0.0 

E-00 

psK 

0.0 

E-00 

QK 

0.0 

E-00 

RK 

Tati 


( Cor'. ' 


c-C) 


9  7 


Flaud  et.  al .  ozone  (668)  010  state  Watson  constan 


684.613410 

E+00 

3.53817681 

E+00 

0.418829269 

E+00 

0.371940160 

E+00 

0.223517 

£-03 

■0.ie5794 

E  -  05 

0.409396 

E  -  06 

0.35057 

E  -  05 

0.602338 

E  -  07 

0.44078 

E-07 

■0.1632 

E  -  08 

■0.92 

E- 12 

0.315 

E- 12 

0.272 

E  -  06 

0.76 

E-il 

0.1500 

E- 12 

0.523 

E-il 

0.95 

E- 13 

0.0 

E-  00 

0.0 

E-  00 

0.0 

E  -  00 

0.0 

L-00 

0.0 

E-00 

c.o 

E  00 

c.o 

E  00 

0.0 

E-00 

c.o 

E-00 

0.0 

E-00 

0.0 

E-00 

0.0 

E-00 

Vibrational  energy 
A 
B 
C 

DLK 

DDK 

DU 

bsK 

dsJ 

HLK 

HLKJ 

HUK 

HU 

hsK 

hsJK 

hsJ 

ALK 

ALKKJ 

AUK 

ALK 

AU7K 

AisK 

AlsJK 

AisKJ 

Al  s  J 

PK 

PKKJ 

psK 

QK 

RK 


Flaud  et.  al .  ozone  (668)  100  state  Watson  constants 


1087.321 

E+00 

Vibrational  energy 

3.492163 

E+00 

A 

0.4177232 

E+00 

B 

0.371583 

E+00 

C 

0.21002 

£-03 

DLK 

0.21110 

E-05 

DUK 

0.41042 

E-06 

DU 

0.367 

E-05 

dsK 

0.5909 

E  -  07 

dsJ 

0.465 

E-07 

HLK 

0.171 

E-08 

HLKJ 

0.300 

E- 11 

HUK 

0.237 

E- 12 

HU 

0.230 

E-08 

hsK 

•0.113 

E  - 10 

hsJK 

C.160 

E- 12 

hsJ 

0.0 

E-00 

ALK 

0.0 

E-00 

ALKKJ 

0.0 

E-00 

AUK 

0.0 

E-00 

ALK 

0.0 

E-00 

AUJK 

0.0 

E-00 

AisK 

0.0 

E-00 

AlsJK 

0.0 

E  00 

AlsKJ 

0.0 

E  00 

AlsJ 

0.0 

E-00 

PK 

o.c 

E-00 

PKKJ 

0.0 

E-00 

psK 

0.0 

E-00 

QK 

0.0 

E-00 

RK 

Table  20.  (Continued) 
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Flaud  et.  ai..  ozone  (668)  001  state  Watson  constan 


1031.144 

E+00 

Vibrational  energv 

3.437126 

E+00 

A 

0.41623161 

E+00 

B 

0.370827 

E+00 

C 

0.20100 

E-03 

dij: 

-0.18600 

E-05 

DUK 

0.41252 

E-  06 

DLJ 

0.2265 

E-05 

asK 

0.6396 

E  -  07 

ds  J 

0.346 

E  -  07 

HLK 

-0.171 

E  -  08 

HLKJ 

-0.300 

E- 11 

HLJK 

0.237 

E- 12 

HLJ 

0.230 

E-08 

hsK 

-0.113 

E- 10 

hsJK 

0.160 

E- 12 

hsJ 

0.0 

E-00 

ALK 

0.0 

E-  00 

ALKKJ 

0.0 

E-00 

AUK 

0.0 

E-00 

ALK 

0.0 

E-00 

AU.'!' 

o 

O 

E-00 

Al  sK 

O 

O 

E-00 

AlsJK 

o 

o 

E  00 

AlsKJ 

0.0 

E  00 

AlsJ 

0.0 

E-00 

PK 

0.0 

E-00 

PKKJ 

0.0 

E-00 

psK 

0.0 

E-00 

QK 

0.0 

E-00 

RK 

Table  20.  (Continued) 
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Isotopic  species  AFGL  notation 


T(K) 

666 

686 

668 

75.0 

423.57301 

442.04036 

452 . 18747 

100.0 

651.93939 

680.37119 

696.01106 

125.0 

911.29072 

951.06373 

972.96670 

150.0 

1199.11158 

1251.55313 

1280.46761 

175.0 

1514.51893 

1581.01883 

1617 . 70104 

200.0 

1857.98773 

1940.05606 

1985.26019 

225.0 

2231.08186 

2330.30476 

2384.75109 

250.0 

2636.19739 

2754.06241 

2818.40078 

275.0 

3076.33884 

3213.93102 

3288.70960 

300.0 

3554.93599 

3712.52993 

3798.17253 

325.0 

4075.69798 

4252.28448 

4349.07628 

350.0 

4642.49715 

4835.28585 

4943.36631 

375.0 

5259.27411 

5463.21039 

5582.57170 

400.0 

5929.95802 

6137.28490 

6267.77471 

*  ‘ible  <!l.  Partition  Suns  from  a  "Complete"  Set  of  0 
Energy  Levels  3 


101 


T(K) 

Q(up)  % 

Err 

Q(actual) 

%  Err 

Q(dowri) 

isotope 

100.0 

653.89971 

.30 

651.93939 

-  .  30 

649.98887 

200.0 

1863.73002 

.  31 

1857.98773 

-  .  31 

1852 .27591 

296.0 

3487.31957 

.  33 

3475.62583 

-  .  34 

3464 . 00203 

300.0 

3566.94909 

.  34 

3554.93599 

-  .  34 

3542 . 99506 

400.0 

5952.39804 

.  38 

5929.95802 

-  .  38 

5907.66662 

isotope 

100.0 

682.41710 

.  30 

680.37119 

-  .  30 

678 . 33551 

200.0 

1946.06082 

.31 

1940.05606 

-  .  31 

1934.08323 

296.0 

3642.21345 

.  33 

3630.04803 

-  .  34 

3617.95402 

300.0 

3725.01776 

.  33 

3712.52993 

-  .  34 

3700.11563 

400.0 

6159.67172 

.36 

6137 . 28490 

-  .  36 

6115.03544 

isotope 

100.0 

698.10417 

.  30 

696.01106 

-  .  30 

693 . 92840 

200.0 

1991.40759 

.31 

1985.26019 

-  .31 

1979.14548 

296.0 

3726.35560 

.33 

3713.92951 

-  .  33 

3701 . 57603 

300.0 

3810.92570 

.33 

3798.17253 

-  .  34 

3785.49411 

400.0 

6290.51425 

.36 

6267.77471 

-  .  36 

6245 . 17363 

T  3.  S*I  G  2  l  . 


Error  Analysis  for  the 


Isotopic  Species  of  0-, 
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Temp(K)  Error  Qactuai 


100. 

.0000 

651.9394 

200. 

- .0008 

1857.9880 

225. 

- .0008 

2231.0820 

250. 

- .0002 

2636.1980 

275. 

.0012 

3076 . 3390 

300. 

.0035 

3554.9360 

325. 

.0066 

4075.6980 

350. 

.0101 

4642.4970 

375  . 

.0135 

5259.2740 

395. 

.0155 

5791.3220 

^tocal 

^vib 

^rot 

651.9395 

1.0000 

651.9119 

11858.0030 

1.0074 

1844.3030 

2231.0990 

1.0136 

2201.1150 

2636.2020 

1.0224 

2578.5150 

3076.3010 

1.0339 

2975. <4^20 

3554.8120 

1.0483 

3391.1030 

4075.4300 

1.0656 

3824.6350 

4642.0260 

1.0858 

4275.3910 

5258.5630 

1.1088 

4742.7660 

5790.4260 

1.1291 

5128.2660 

Tat: 


-  c  s t  of  Approximation  Q  . 

'actual 

the  Principle  Species  of  0^ 


lvib  xrot 


for 


103 


1 

Toth's  Effective  rotational 

constants 

for  the 

0000  state  of  446  N20 

1  1 

91  120 

0.0000000000 

E  00 

G 

0.4190110012 

E  00 

B 

1.76091933 

E-07 

D 

-0.1652924 

E-13 

H 

0.00000 

E-18 

L 

1  Toth's  Effective  rotational 

constants 

for  the 

0110E 

state  of 

446  N20 

1  2E 

80  120 

588.768 

E  00 

G 

0.4191779259 

E  00 

B 

1.78324547 

E-07 

D 

-0.17146 

E-13 

H 

0.00000 

E-18 

L 

I  Toth's  Effective  rotational 

constants 

for  the 

0110F 

state  of 

446  N20 

1  2F 

80  120 

588.768 

E  00 

G 

0.4199698451 

E  00 

B 

1.79303055 

E-07 

D 

-0.17668 

E-13 

H 

0.00000 

E-18 

L 

1  Toth's  Effective  rotational 

constants 

for  the 

0200  state  of  446  N20  | 

1  3 

78  120 

1168 . 132-3 

E  00 

G 

0.419920957 

E  00 

B 

2.49403822 

E-07 

D 

30.82612 

E-13 

H 

-18.9735 

E-18 

L 

|  Toth's  Effective  rotational 

constants 

for  the 

0220E 

state  of 

446  H20 

1  4E 

78  120 

1177.7446 

E  00 

G 

0.420125256 

E  00 

B 

1.196793 

E-07 

D 

-29.5680 

E-13 

H 

-24.48239 

E-18 

L 

R  Toth's  Effective  rotational 

constants 

for  the 

02  2  OF 

state  of 

446  N20 

1  4F 

78  120 

1177.7446 

E  00 

G 

0.420126260 

E  00 

B 

1.8180 

E-07 

D 

0.9510 

E-13 

H 

0.00000 

E-18 

L 

Table  2  u  .  V 

ilrat ional-Fc 

ta  l  2''  nal 

Han i It 

o  r ;  i  t  r. 

Ccr.star 

ts 

or  from  T 

L. 

oti 

10‘* 

Toth's  Effective  rotational  constants  for  the  1000  state  of  446  N20 


1  5 

1284.9 
0.417255097 
1.7261943 
1.263279 
3.492826 


76  120 

E  00 
E  00 
E-07 
E-13 
E-18 


Toth's  Effective  rotational 
1  6E  68  120 

1749.0651  E  00 

0.419583988  E  00 

2.109928  E-07 

11.75046  E-13 

-12.20395  E-18 

Toth's  Effective  rotational 
1  6F  68  120 


1749.0651 

0.421079177 

2.178319 

-3.200163 

14.50675 


E  00 
E  00 
E-07 
E-13 
E-18 


Toth's  Effective  rotational 
1  7E  67  120 

1766.912  E  00 

0.420667057  E  00 

1.593325  E-07 

-29.1628  E-13 

0.0  E-18 

Toth's  Effective  rotational 
1  7F  67  120 

1766.912  E  00 

0.4206713799  E  00 
1.6806495  E-07 

29.5664  E-13 

0.0  E-18 

Toth's  Effective  rotational 


G 

B 

D 

H 

L 

constants 

G 

B 

D 

H 

L 

constants 

G 

B 

D 

H 

L 

constants 

G 

B 

D 

H 

L 

constants 

G 

B 

D 

H 

L 

constants 


for  the  0310E  state  of  446  N20 


for  the  031 OF  state  of  446  N20 


for  the  0330E  state  of  446  N20 


for  the  03 3 OF  state  of  446  N20 


for  the  0330E  state  of  446  N20 


1  7E 
1766.912 
0.420667057 
1.593325 
-29.1628 
0.0 


1  2 
E  00 
E  00 
E-07 
E-13 
E-18 


( Continued ) 


105 


Toth's  Effective  rotational  constants  for  the  03 3 OF  state  of  446  N20 
1  7F  1  2 

1766.912  E  00  G 

0.4206713799  E  00  B 

1.6806495  E-07  D 

29.5664  E-13  H 

0.0  E-18  L 

Toth's  Effective  rotational  constants  for  the  111CE  state  of  4<^ 

1  8E  66  120 

1880.2657  E  00  G 

0.4174647588  E  00  B 

1.74899344  E-07  D 

1.138277  E-13  H 

0.00000  E-18  L 

Toth's  Effective  rotational  constants  for  the  1110F  sta.  -  cf  -14>- 


1  8F 

66  120 

1880.2657 

E  00 

G 

0.4183730749 

E  00 

B 

1.72005845 

E-07 

D 

2.246281 

E-13 

H 

0.00000 

E-18 

L 

10  6 


ur-:) 


1 


Toth's  Effective  rotational  constants  for  the  0000  state  of  456  N20 
2  1  78  120 

0.0000000000  E  00  G 

0.4189818102  E  00  B 

1.76326469  E-07  D 

0.0  E-13  H 

0.00000  E-18  L 

Toth's  Effective  rotational  constants  for  the  0110E  state  of  456  N20 
2  2E  68  120 

588.768  E  00  G 

0.4191089166  E  00  B 

1.78583283  E-07  D 

0.0  E-13  H 

0.00000  E-18  L 

Toth's  Effective  rotational  constants  for  the  0110F  state  of  456  N20 


2  2F 
588.768 
0.4199186417 
1.79446290 
0.0 

0.00000 


68  120 
E  00 
E  00 
E-07 
E-13 
E-18 


G 

B 

D 

H 

L 


Toth's  Effective  rotational  constants  for  the  0200  state  of  456  N20 


2  3 

1  120 

1168.1323 

E  00 

G 

0.4199204111 

E  00 

B 

2.8246994 

E-07 

D 

58.90136 

E-13 

H 

0.0 

E-18 

L 

Toth's  Effective  rotational 

cons 

2  4E 

1  120 

1177.7446 

E  00 

G 

0.42003233 

E  00 

B 

0.9542 

E-07 

D 

0.0 

E-13 

H 

0.0 

E-18 

L 

for  the  0220E  state  of  456  N20 


Table  24 . 


(Continued ) 
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Toth's  Effective  rotational  constants  for  the  0220F  state  of  456  N20 
2  4F  1  120 


1177.7446  E  00  G 

0.42003483  E  00  B 

1.8131  E-07  D 

0.0  E-13  H 

0.00000  E-18  L 


Toth's  Effective  rotational  constants  for  the  1000  state  of  456  N20 


25  1  120 


1284.9  E  00  G 

0.417127383  E  00  B 

1.75528298  E-07  D 

0.60935  E-13  H 

0.0  E-18  L 


Toth's  Effective  rotational  constants  for  the  1110E  state  of  456  N20 
2  8E  1  120 


1880.2657  E  00  G 

0.4172733118  E  00  B 

1.77554100  E-07  D 

0.562150  E-13  H 

0.00000  E-18  L 


Toth's  Effective  rotational  constants  for  the  1110F  state  of  456  N20 


2  8F  1  120 


1880.2657  E  00  G 

0.4181654622  E  00  B 

1.76635130  E-07  D 

2.288250  E-13  H 

0.00000  E-18  L 


Table  2-.  (Continued) 
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Toth's  Effective  rotational 
3  1  79  120 


constants  for  the  0000  state  of  546  N20 


0.0000000000  E  00 

0.404857965155  E  00 

1.64294216  E-07 

1.867536  E-13 

0.00000  E-18 

Toth's  Effective  rotational 


3  2E 
588.768 
0.4050372717 
1.65680840 
0.0 

0.00000 


69 


120 
E  00 
E  00 
E-07 
E-13 
E-18 


Toth's  Effective  rotational 
3  2F  69  120 

588.768  E  00 

0.40578110957  E  00 

1.66742291  E-07 

0.0  E-13 

0.00000  E-18 

Toth's  Effective  rotational 
33  1  120 

1168.1323  E  00 

0.40571242618  E  00 

2.224294126  E-07 

21.26045  E-13 

0.0  E-18 

Toth's  Effective  rotational 


3  4E 
1177.7446 
0.405951230 
1.21815 
0.0 
0.0 


1  120 
E  00 
E  00 
E-07 
E-13 
E-18 


Toth's  Effective  rotational 


3  4F 
1177.7446 
0.405950430 
1.69110 
0.0 

0.00000 


1  120 
E  00 
E  00 
E-07 
E-13 
E-18 


G 

B 

D 

H 

L 

constants 

G 

B 

D 

H 

L 

constants 

G 

B 

D 

H 

L 

constants 

G 

B 

D 

H 

L 

constants 

G 

B 

D 

H 

L 

constants 

G 

B 

D 

H 

L 


for  the  0110E  state  of  546  N20 


for  the  0110F  state  of  546  N20 


for  the  0200  state  of  546  N20 


for  the  0220E  state  of  546  N20 


for  the  02 2 OF  state  of  546  N20 


Table 
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Toth's  Effective  rotational  constants  for  the  1000  state  of  546  N20 
35  1  120 

1284.9  E  00  G 

0.403263631305  E  00  B 

1.5980072  E-07  D 

3.702650  E-13  H 

0.0  E-18  L 

Toth's  Effective  rotational  constants  for  the  0310E  state  of  546  N20 


3  6E 
1749.0651 
0.405410701 
1.9139457 
0.0 
0.0 


L  120 
E  00 
E  00 
E-07 
E-13 
E-18 


G 

B 

D 

H 

L 


Toth's  Effective  rotational  constants  for  the  0310F  state  of  546  N20 
3  6F  1  120 

1749.0651  E  00  G 

0.4067990584  E  00  B 

1.99573470  E-07  D 

0.0  E-13  H 

0.0  E-18  L 

Toth's  Effective  rotational  constants  for  the  1110E  state  of  546  N20 
3  8E  1  120 

1880.2657  E  00  G 

0.40348977155  E  00  B 

1.61182100  E-07  D 

0.657200  E-13  H 

0.00000  E-18  L 

Toth's  Effective  rotational  constants  for  the  1110F  state  of  546  N20 
3  8F  1  120 

1880.2657  E  00  G 

0.40435854185  E  00  B 

1.57574110  E-07  D 

0.244250  E-13  H 

0.00000  E-18  L 


fable  21-.  CCcn~i:,ue!) 
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Toth's  Effective  rotational  constants  for  the  0000  state  of  448  N20 


4  1  78  120 

0.0000000000  E  00 

0.39557789529  E  00 

1.58345694  E-07 

0.0  E-13 

0.00000  E-18 

Toth's  Effective  rotational 
4  2E  68  120 

588.768  E  00 

0.39575970106  E  00 

1.60023584  E-07 

0.0  E-13 

0.00000  E-18 

Toth's  Effective  rotational 


4  2F  68  120 


588.768  E  00 

0.396471172725  E  00 
1.61077378  E-07 

0.0  E-13 

0.00000  E-18 


Toth's  Effective  rotational 


43  1  120 

1168.1323  E  00 

0.396278094345  E  00 
2.05540250  E-07 

11.96215  E-13 

0.0  E-18 

Toth's  Effective  rotational 
45  1  120 

1284.9  E  00 

0.394111321136  E  00 
1.49508691  E-07 

3.498220  E-13 

0.0  E-18 


G 

B 

D 

H 

L 

constants  for  the  0110E  state  of  448  N20 

G 

B 

D 

H 

L 

constants  for  the  0110F  state  of  448  N20 

G 

B 

D 

H 

L 

constants  for  the  0200  state  of  448  H20 

G 

B 

D 

H 

L 

constants  for  the  1000  state  of  448  N20 

G 

B 

D 

H 

L 


Tab: 
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Toth's  Effective  rotational  constants  for  the  0310E  state  of  448  N20 


4  6E 
1749.0651 
0.3959890999 
1.836545 
0.0 
0.0 


1  120 
E  00 
E  00 
E-07 
E-13 
E-18 


Toth's  Effective  rotational  constants  for  the  0310F  state  of  448  N20 


4  6F 
1749.0651 
0.3972774667 
1.93453325 
0.0 
0.0 


1  120 
E  00 
E  00 
E-07 
E-13 
E-18 


Toth's  Effective  rotational  constants  for  the  1110E  state  of  448  N20 


4  8E 
1880.2657 
0.3943490062 
1.52755220 
0.0 
0.0 


1  120 
E  00 
E  00 
E-07 
E-13 
E-18 


Toth's  Effective  rotational  constants  for  the  1110F  state  of  448  K20 


4  8F 
1880.2657 
0.3952190532 
1.47157125 
0.0 
0.0 


1  120 
E  00 
E  00 
E-07 
E-13 
E-18 


Toth's  Effective  rotational  constants  for  the  0000  state  of  447  N20 
5  1  78  120 


0.0000000000  E  00  G 

0.406672154476  E  00  B 

1.66397195  E-07  D 

0.0  E-13  H 

0.0  E-18  L 
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Toth's  Effective  rotational  constants  for  the  0110E  state  of  447  N20 


5  2E 

68  120 

588.768 

E  00 

G 

0.4068471 

E  00 

B 

1.6836 

E-07 

D 

0.0 

E-13 

H 

0.0 

E-18 

L 

Toth's  Effective  rotational 

constants 

5  2F 

68  120 

588.768 

E  00 

G 

0.4075991 

E  00 

B 

1.7139 

E-07 

D 

0.0 

E-13 

H 

0.0 

E-18 

L 

Toth's  Effective  rotational 

constants 

5  3 

1  120 

1168.1323 

E  00 

G 

0.40748554120 

E  00 

B 

2.222921200 

E-07 

D 

13.82892 

E-13 

H 

0.0 

E-18 

L 

Toth's  Effective  rotational 

constants 

5  5 

1  120 

1284.9 

E  00 

G 

0.4050539924 

E  00 

B 

1.6100012 

E-07 

D 

3.37195 

E-13 

H 

0.0 

E-18 

L 

Tatle 


( Cor.*'  inu^d  ) 
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ivlow 


vibrational  state 


ivlow 


vibrational  state 


1  0000  5  1000 

2  0110E,  001 10F  6  0310E,  0310F 

3  0200  7  0300E,  0300F 

4  0200E,  0200F  8  1110E,  1 1 1  OF 

For  N2O  the  energy  levels  can  be  expressed  by  the  formula: 

Evr=  G  +  B  -  DX2+  HX3+  LX4  ^5) 

where  X  =  J(J  +  1)  and  G,  B,  D,  H  and  L  are  effective  constants  given  in 
reference  41  (see  also  Table  24).  The  missing  J  levels  for  each 
isotopic  species  and  each  vibrational  state  are  given  below  in  Table 
25  while  the  band  origins  for  the  vibrational  states  from  Toth42  are 
given  in  Table  26. 

Table  25.  Last  J-Level  for  N2O  Data  From  the  AFGL  Database 

Vibrational  State 


isotopic 

species 

0000 

0110 

0200 

0220 

1000 

0310 

0300  1110 

446 

91 

80 

78 

78 

76 

68 

67  66 

456 

78 

68 

546 

79 

69 

448 

78 

68 

447 

77 

60 

lit- 


Table  26.  Band  Origins  for  the  Various  Vibrational  Levels 


Vib. 

Level 

Band  Origin 
for  Isotopic 

in  cm*1 
Species 

446 

456 

546 

448 

447 

0100 

588.768 

575.5 

585.320 

584.1 

586.3 

0200 

1  168.1  323 

1  144.3334 

1  159.9717 

1  155.1  399 

0220 

1  177.7446 

1000 

1284.9 

1269.8920 

1246.8846 

1264.7047 

03  10 

1749.065  1 

0300 

1766.912 

1110 

1880.2657 

1860.1913 

1862.7670  1839.9361 

Some  4000  energy  levels  were  calculated  to  fill  the  N2O  energy 
array.  Because  of  the  time  necessary  for  file  transfer  to  and  from  the 
VAX  an  energy  cutoff  of  5000  cm*1  was  enforced  bringing  the 
number  of  calculated  levels  to  2888.  The  complete  set  of  energy 
levels  (to  5000  cm*1)  was  then  used  to  calculate  the  partition  sums 
for  N2O.  In  the  error  analysis  part  of  the  calculations,  the 
recalculated  (high  J)  levels  are  flagged  by  a  negative  sign  and  a  shift 
of  0.5%  is  used.  For  all  other  levels  the  0.2%  shift  was  used.  Table  27 
lists  the  partition  sums  for  N2O.  The  error  analysis  gave 
approximately  a  0.5%  error  at  500  K  for  the  546  species  as  the 
largest  error. 


5.4.1  CO 


The  total  internal  partition  sums  for  CO  were  calculated; 
however,  the  polynomial  fit  has  not  yet  been  done.  For  this  molecule 


T(K)  Isotopic  species  AFGL  notation 


446 

456 

546 

448 

47  7 

75.0 

124.77949 

124.78836 

129.13001 

132.15226 

128  .  548  3  3 

100.0 

166.33257 

166 . 35749 

172.13874 

176.17161 

171 . 3624 

200.0 

342.09902 

343.05349 

354.26476 

362 . 50801 

352.4321 7 

225.0 

391.78094 

393.28109 

405.81071 

415.03689 

403.4433 3 

250.0 

445.02930 

447 . 16813 

461 .06764 

471.07548 

457 . 807  34 

270.0 

490.54599 

493.22996 

508.29755 

518.68522 

503.95600 

296.0 

554.01800 

557  .  39564 

574.13339 

584.53580 

567 . 70333 

300.0 

564.24195 

567 .71963 

584 .73393 

595.07935 

577.89900 

325.0 

631.08477 

635.10712 

654.00040 

663 . 55749 

644.02468 

350.0 

703.26770 

707.61797 

728 . 70833 

736.61184 

714.35619 

400.0 

864.98315 

868.84845 

895.63565 

897.03323 

867.81020 

500.0 

1265.20127 

1259 . 64420 

1305.56265 

1278.03489 

1225.49795 

Table  27. 


Partition 


nr. 


cf  Energy  Levels; 


the  energy  levels  were  calculated  using  the  Dunham  expansion.  The 
original  program  and  Dunham  coefficients  were  supplied  to  us  by  R. 
Tipping,  University  of  Alabama.  Because  the  number  of  levels  is 
relatively  small  we  have  calculated  all  levels  from  J  =  0  to  74  for  the 
ten  lowest  vibrational  levels  (v  =  0,  9).  This  generates  energies  to  = 
28,000  cm'l.  The  calculations  were  done  for  the  four  isotopic  species 
of  CO  on  the  database. 


5.5  Methane 


For  the  methane  molecule,  complete  energy  level  listings 
for  the  ground  states  of  the  and  species  were  obtained 

from  Dr.  Linda  Brown,  Jet  Propulsion  Laboratory.  The  energies  are 
given  for  J  up  to  31  (~5000.  cnr*).  These  were  merged  with  the 
upper  vibrational  state  energies  from  the  database  to  get  a 
"complete"  set  of  energies.  The  partition  sums  were  calculated  for 
these  species  and  are  given  in  Table  28.  In  Figure  21  a  plot  of  the 
partition  sum  for  each  species  of  carbon  vs.  energy  at  400  K  is 
presented.  We  note  the  sums  have  converged.  The  polynomial 
coefficients  for  the  temperature  expression  were  then  calculated  for 
this  molecule. 


5.6  Nitrogen  Dioxide 

For  nitrogen  dioxide  only  the  principal  isotopic  species  is 
present  on  the  database,  14N1602-  The  coupling  of  the  rotational 
wave  functions  and  the  nuclear  statistics  for  the  symmetric  oxygen 
atoms  for  this  molecule  is  the  same  as  for  symmetric  O3,  i.e.  the  only 
allowed  levels  are  for  Ka  +  Kc  +  V3  even.  All  energy  levels  were 
determined  via  theory.  The  molecular  constants  were  obtained  for 
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L-  CO 


211  isotope 
311  isotope 

:  i  vir-;  .  . .  Partition  Sun  vs.  Energy  for  ""CH  and  A  'CH 


the  first  eleven  vibrational  states  of  N02-  The  constants  were  taken 
from  reference  43  for  the  000,  010,  020,  100,  001,  030  and  011 
bands  and  from  reference  44  for  the  120,  101,  130  and  the  111 

bands.  All  energies  from  J  =  0  to  75  were  calculated  and  those  less 
than  5000  cnr*  were  retained. 

Nitrogen  dioxide  has  an  unpaired  electron  which 
complicates  its  structure.  The  rotational  quantum  numbers  are 

labeled  by  N,  Ka  and  Kc  (note  Ka  and  Kc  can  be  combined  into  one 
quantum  number  x  =  Ka-Kc)  where  N  is  the  angular  momentum, 
rotational  and  electron  orbital,  and  Ka  and  Kc  are  the  prolate  and 
oblate  limit  of  K.  The  electron  spin  (S  =  1/2)  and  rotational  angular 
momentum  (N)  couple  to  give  the  total  electron  and  rotational 

angular  momentum  J  =  N  ±S  =  N±  1/2,  thus  the  spin-rotation 

interaction  splits  the  energy  levels  into  two  components.  The  energy 
can  be  expressed  as  a  rotational  energy,  Er  given  by  diagonalizing  the 
Watson  Hamiltonian,  and  a  contribution  due  to  the  spin-rotation 
interaction.  The  spin-rotation  interaction  terms  determined  up  to 
second  order  are: 


(eaa-e)Ka  +£N(N-H)  + y-(C~- Ecc)  5,  K  N(N+1) 
1  4  '  ' 


he  2(N+1) 


+  n  ^aa~E^  k2  1 

+  T'aaaa  —  *^a  *  . 

4  R  2 

(N+l) 


for  J  =  N  +  1/2,  and 


, -  (eaa-E)K l  -f  e  N(N+1)  +  y  (-b-  1C^-S1  N(N+1) 

hr  |  dJ  4 


he  2N| 


<a  aaa  a 


Ka  -  (£at,-  kJ  1  -  — 
4B  N2 


for  J  =  N  -  1/2  (N>0) 

In  these  equations  e  =  (1/2)  (£bb'ecc)i  ®  =  0/2)(B+C);  y  =  -1  for  N  odd 
and  y=  +1  for  N  even.  The  spin-rotation  constants  are  given  in  terms 
of  the  eaa,  £bb>  ecc  ancl  ’laaaa-  From  the  papers  listed  above,  spin- 
rotation  constants  are  available  for  the  following  bands:  000,  010, 
120,  101,  130,  and  111.  For  the  other  bands,  020,  100,  001,  030  and 
011,  only  the  rotation  energy  was  calculated  and  a  splitting  factor  of 
2  was  multiplied  times  these  levels  in  the  partition  sums,  i.e.  the 
degeneracy  factor. 

We  started  by  calculating  the  rotational  contribution  to 
the  energy  levels.  All  states  for  N  =  0  to  75  were  considered  with  the 
cutoff  at  5000  cm*1.  The  calculations  had  some  difficulties:  for 

lower  values  of  N  the  calculations  matched  the  published  values, 
however,  extrapolation  to  higher  values  was  clearly  incorrect.  For 
example,  the  state  44  6  38  is  present  on  the  database  and  has  an 
energy  of  1105.6  cm'1,  the  calculated  value  from  diagonalizing  the 

Watson  Hamiltonian  was  402.9  cm'1.  From  the  calculations  there 
were  many  other  examples  of  levels  that  were  incorrect,  some  wild 
extrapolations,  negative  energies,  etc.  In  order  to  calculate  the 

partition  sums  accurate  to  one  percent,  these  levels  must  be 

evaluated.  We  looked  at  several  ways  to  estimate  the  missing  levels, 
e.g.  filling  in  with  levels  that  were  present  on  the  database  and 


trying  fits  for  the  other  missing  levels.  For  NCH,  the  rotational 
constants  B  and  C  are  approximately  equal  to  give  a  near  prolate 
symmetric  top  but  the  limiting  energy  formula  does  not  give 
accurate  enough  results  for  the  calculations.  We  also  investigated  the 
change  in  energy  with  respect  to  Ka  and  again  this  was  of  little  help 
for  prediction  of  the  energies  that  were  needed. 

We  then  considered  the  fitting  procedure  that  is  normally 
used  to  give  the  Watson  Hamiltonian  constants.  As  one  goes  up  in 
the  terms  of  the  Hamiltonian  the  constants  get  smaller  and  smaller. 
However,  these  multiply  angular  momentum  to  larger  and  larger 
powers.  Thus,  if  one  of  the  upper  constants  is  off  by  a  small  amount, 
the  calculations  diverge  as  J  increases.  We  then  took  the  ground 
state  Watson  constants  from  Camy-Peyret  et.  al.4^  (this  publication 
references  the  work  of  W.  Lafferty,  N.B.S.  for  the  ground  state 
Hamiltonian  constants)  and  set  the  three  highest  constants,  %N- 
and  Lk,  to  zero  (these  are  angular  momentum  to  the  sixth,  sixth,  and 
eighth  powers,  respectively).  The  calculations  were  repeated  with 
these  constants  and  compared  with  the  published  values  (those  on 
the  HITRAN  database).  Table  29  gives  examples  for  some  levels  that 
were  in  error  before. 

Table  29.  Comparison  of  Calculated  Energies  in  cm"'  with 

Published  Values 


State 

E(calc) 

E(pub) 

%  diff 

29  10  20 

1099.8 

1 102. 

0.2 

12  12  0 

1102.1 

1 109. 

0.62 

44  6  38 

1 105.6 

1105.6 

0.0 

51  151 

1 106.8 

1 106.7 

-0.005 

The  agreement  is  very  good  and  it  appears  that  this  was  the 
problem  We  tested  eliminating  only  the  highest  constant  Lk  (this  is 
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the  Pz8  term)  and  calculating  the  energy.  The  results  agree  very  well 
with  the  published  values  so  we  feel  it  was  only  the  Lk  constant  that 
was  responsible  for  the  error  in  the  diagonalizations.  This  constant 
was  removed  from  all  the  sets  of  Hamiltonian  constants  and  the 
energies  were  calculated  for  all  levels  of  NO2;  this  generated  15,984 
rotational  energy  levels.  In  Table  30,  the  Watson  Hamiltonian 
constants  used  in  the  calculations  are  given.  These  are  from  the 
papers  cited  above  with  the  modifications  discussed  in  the  text.  We 
wrote  a  program  to  calculate  the  spin-rotation  splitting  to  the  energy 
levels  for  the  six  bands  that  constants  are  available  for.  These 
constants  are  given  in  Table  31. 

Table  31.  Spin-Rotation  Interaction  Constants  for  NO2 


Band 

e  =  l/2(ebb-ecc) 

£aa-£ 

£bb-Ecc 

Haaaa 

000 

-0.14601935  x  10  2 

0.18180028 

0.34220076x10-2 

-0.1710183x10-3 

010 

-0.15193  x  10*2 

0.199308 

0.384x10-2 

-0.20074x10-3 

120 

-0.16733  x  10-2 

0.21940 

0.446x10-2 

-0.24054x10-3 

101 

-0.15513  x  10-2 

0.17648 

0.3708x10-2 

-0.1733x10-3 

130 

-0.17321  x  10-2 

0.236906 

0.48675x10-2 

-0.26952x10-3 

111 

-0.16189  x  10-2 

0.193993 

0.41 165x10*2 

-0.20228x10-3 

Applying  this,  a  total  of  24,384  energy  levels  were  generated  for 
NO2,  16804  states  which  have  the  spin-rotation  interaction  included 
and  7580  states  which  do  not  and  have  an  additional  factor  of  2  in 
the  partition  sum  degeneracy  factor.  The  non-symmetrical  nitrogen 
nucleus  I(14N)  gives  a  nuclear  statistical  factor  of  gj  =  2(1  +  1)  =  3. 
The  partition  sums  were  calculated  over  the  temperature  grid  to 
allow  fitting  the  polynomial  to  Q(T).  The  polynomial  fit  reproduced 
the  sums  to  better  than  0.4%.  Comparison  of  our  calculated  value  at 
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225K  with  the  value  given  on  the  JPL  Catalogue  gives  Q(ULCAR)  = 
8835  and  Q(JPL)  =  8750. 

5.7  Approximate  Partition  Sums 

As  part  of  the  work  on  calculating  the  total  internal 
partition  sums  for  the  molecules  on  the  database,  we  tested  an 
approximation  that  is  often  used  for  evaluating  the  partition  sum. 
The  approximation  is  to  replace  the  total  internal  partition  sum  by  a 
product  of  the  vibrational  partition  sum  and  the  rotational  partition 
sum,  i.e.  Q(T)  =  Qv(T)*Qr(T).  Rather  than  doing  a  direct  sum  on 
twenty  to  thirty  thousand  energy  levels,  two  smaller  direct  sums  are 
done.  The  vibrational  sum  is  usually  over  5  to  10  terms  and  the 
rotational  sum  can  be  over  several  hundred  to  a  thousand  or  so 
levels  or  can  be  approximated  by  analytical  formula  or  by  classical 
values. 


The  approximation  made  would  be  strictly  true  if  the 
shape  of  each  upper  vibrational  potential  curve  was  the  same  as  the 
ground  state  potential  curve.  This  is  not  the  case,  however,  so  we 
must  ask  the  question  "How  much  does  the  change  in  the  shape  of 
the  potential  curve  affect  the  rotational  energies  and  does  this  have  a 
noticeable  effect  on  the  partition  sums?"  What  we  ust  do  to 
answer  this  is  to  consider  the  ro-vibrational  energies  since  the 
assumption  comes  from  writing  EVr  =  Ev  +  Er.  For  most  of  the  levels 
at  the  bottom  of  the  potential,  the  difference  in  energy  causes  little 
change  in  the  calculated  partition  sum  (<0.1%)  over  temperatures 
below  250-300  K.  Above  these  tempeiatures,  upper  levels  become 
important  and  we  expect  to  see  some  deviation  between  the  two 

partition  sums.  We  desire  to  know  how  great  this  deviation  can  be. 

This  will  be  important  for  the  larger  systems  where  molecular 

constants  are  not  available  for  many  of  the  low  lying  vibrational 
bands.  It  is  for  these  systems  where  the  approximation  must  be 
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made.  This  study  will  help  determine  the  accuracy  of  the  partitions 
sums  calculated  using  Qv*Qr- 

We  have  tested  the  approximation  for  H2O,  CO2  and  O3. 
For  water  vapor  the  calculations  were  done  for  the  principal  isotopic 
species  and  for  monodeuterated  water.  The  direct  summation  for  the 
vibrational  partition  sums  was  over  8  and  10  vibrational  levels 
respectively.  The  values  of  the  vibrational  energy  levels  are  given  in 
Table  32  for  all  the  species  we  have  tested  except  CO2  for  which  all 
of  the  values  published  were  used  (see  L.S.  Rothman,  Appl. 
Opt. 25, 1795  (1986)).  For  H2O.  CO2  and  O3  the  error  was  of  order 
0.1%  in  the  temperature  range  of  interest.  For  HDO  tne  error  was  low 
until  around  275-300  K  where  the  error  went  to  _1%.  This  is  shown 
in  Figure  22.  Also  given  in  the  figure  is  the  approximation  using  the 
classical  value  for  the  rotational  partition  sum.  Note  this  is  clearly 
greater  than  the  1%  level  we  have  set  for  the  calculations. 

We  decided  to  also  test  the  approximation  that  is  used  in 
the  FASCODE  program21  for  the  partition  sums  as  a  function  of 
temperatures.  In  the  program  the  partition  sum  at  temperature  T  is 
approximated  by  the  product  Qv(T)*Qr(T),  where  QV(T)  is  calculated 
by  direct  summation  and  Qr(T)  is  approximated  by 

Qr(296K)*(T/296K) 1 -5.  In  Figure  23  the  calculated  Q(T)  is  plotted  vs. 
the  FASCODE  approximation  for  HDO.  The  %  error  plot  reveals  that 
the  approximation  starts  around  3%  low  at  70K  and  increases  to 
roughly  0.5%  high  at  420  K.  The  approximation  is  not  too  severe  in 
the  temperatures  of  atmospheric  interest,  this  is,  however,  only  for 
one  species  and  larger  deviations  may  be  observed  for  other 

molecules. 


1  3  >■ 


Temperature  K 


Table  32.  Vibrational  Hnergies  in  cm 


level 

IbO 

10X7 

16C >3 

1 

().() 

0.0 

0.0 

-) 

1  394.7 

1403.4 

700.9 

3 

3  13  1.6 

2723.7 

1042. 1 

4 

3637.1 

2782.0 

1  103. 1 

3 

3  733. 9 

3707.4 

1  399.3 

6 

4  666.8 

4100. 1 

2063.2 

7 

3  233.0 

4  1  43.6 

21  10.8 

X 

3 331. 3 

3363.6 

2193.9 

9 

6413.6 

1796.3 

1  0 

6432.1 

1726.3 

3.8  Nuclear  Statistical  Factors  for  Non-Svmmctric  Nuclei 

We  decided  to  calculate  the  internal  partition  sum 
explicitly.  This  generates  nuclear  statistical  factors  that  cancel  in 
usual  applications  of  the  partition  sums  and  hence  these  terms  arc 
normally  omitted.  It  is  more  for  the  sake  of  clarity  that  they  are 
included  here,  i.e.  rather  than  trying  to  explain  what  factors  we  are 
leaving  out  from  the  sum  for  each  isotopic  variant  we  are  simply 
calculating  the  partition  sum  correctly.  In  order  to  use  these  sums, 
e.g.  calculating  l<  from  S,  these  factors  must  be  used.  We  review 
them  here  for  the  species  we  have  considered  to  date.  The  fact  of 
having  non  symmetric  nuclei  means  that  the  spin  of  the  nucleus  has 
not  been  accounted  for  in  the  coupling  between  rotational 
wavcfunctions  and  the  spins  ot  symmetric  nuclei  (the  usual  gj  terms. 


e.g.  0  or  1  for  asymmetric/symmetric  ozone).  The  nuclear  statistical 
factor  generates  additional  multiplying  terms  of  gj  =  (21+1)  for  each 
non-symmetric  nuclei,  where  I  is  the  spin  of  the  nucleus  in  question. 
We  write  the  internal  partition  sum  as 


Q(T)  =  gj  X  §i  e "  /kT 


were  gj  can  be  taken  out  of  the  sum  because  it  is  the  same  for  all 
states.  In  the  discussion  below  use  is  made  of  the  AFGL  isotopic 
labeling  scheme,  e.g.  161  represents  the  H2l^O  species  of  water 
vapor. 


1 6 1  The  hydrogen  nuclei  are  symmetric  and  couple  with  the 

rotational  wave  functions  to  give  the  3:1  (oddteven 
states)  ratio.  The  nucleus  has  spin  zero  so  the 

additional  factor  is  gj  =  [2*(0)+l]  =  1. 

181  This  is  similar  to  161,  because  has  spin  zero  the  factor 
is  the  same  as  for  161,  gj  =  [2*(0)+l]  =  1. 

1 7 1  Again  the  symmetric  hydrogen  nuclei  couple  with  the 
rotational  wavefunctions  to  give  the  3:1  ratio  (odd:even). 
The  170  nuclei  has  a  spin  of  5/2  to  give  a  factor  of  gj  = 
[2*(5/2)+l]  =  6. 

162  The  hydrogen  nuclei  are  not  symmetric  so  there  is  no 

coupling  with  the  rotational  wavefunctions  (gj  =  1).  The 

nuclear  spins  arc  1(H)  =  1/2,  1(D)  =  1  and  I(!6o)=0  to  give 
a  factor  of  gj  =  12*0/2 )+ 11  [2*(1)+1)  |2*(0)+1 1  =  2*3*  1  =  6. 


CCb 


62  6  There  is  coupling  of  the  spin  zero  oxygens  with  the 

rotational  wavefunctions,  has  zero  spin  therefore  one 
has  gj  =  1. 

63  6  Same  argument  for  the  oxygens  as  above,  has  spin 

1/2  to  give  gj  =  [2*(l/2)+l]  =  2. 

628  No  symmetric  nuclei  so  we  get  n[2Ij+l]  but  all  the  spins 
are  zero  so  gj  =  1. 

627  gj  =  [2*(0)+l ]  [2*(0)+l ]  [2*(5/2)+l]  =  1*1*6  =  6 

63  8  gj  =  [2*(0)+l  1  l2*(l/2)+l]  [2*(0)+l]  =  1*2*1  =  2. 

637  gj  =  (2*(0)+l  1  [2*(l/2)+l]  [2*(5/2)+l ]  =  1*2*6  =  12. 

628  gj  =  [2*(0)+l ]  [2*(0)+l]  [2*(0)+l]  =  1*1*1  =  1. 

7  28  gj  =  2*(5/2)+l I  [2*(0)+l j  [2*(0)+l]  =  6*1*1  =  6. 


666  Symmetric  end  oxygens  couple  with  rotational 
wavefunctions  to  give  the  1:0  (even:odd)  ratio.  The 
central  oxygen  (non-symmetric)  has  spin  zero  to  give  gj  = 
[2*(0)+l ]  =  1. 

66  8  No  coupling  between  rotational  wavefunctions  and  spin 
so  all  three  oxygens  contribute  to  give  gj  =  [2*(0)+l] 
12*(0)+1]  (2*(0)+l  J  =  1*1*1  =  1. 
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68  6  Coupling  for  the  1:0  ratio  as  for  666,  the  central  has 

zero  spin  to  give  gj  =  [2*(0)+l]  =  1. 

n2o 

44  6  There  are  no  symmetric  nuclei  for  this  molecule  thus  one 
has  gj  =  Fl[2Ij+l].  One  has  for  446,  gj  =  [2*(1)+1  ]  [2*(1)+1] 
[2*(0)+l]  =  3*3*1  =  9. 

456  gj  =  [2*(1)+1]  |2*(l/2)+l]  [2*(0)+l]  =  3*2*1  =  6. 

546  gj  =  [2*(l/2)+l]  [2*(1  )+l  ]  (2*(0)+l  1  =  2*3*1  =  6. 

448  gj  =  [2*(1  )+l ]  (2*(1)+H  12*(0)+1]  =  3*3*1  =  9. 

447  gj  =  [2*(1  )+l  ]  [2*(1)+1]  [2*(5/2)+l  ]  =  3*3*6  =  54. 

CO 

26  gj  =  [2*(0)+l]  [2*(0)+lj  =  1*1  =  1. 

3  6  gj  =  [2*(l/2)+l]  [2*(0)+l  ]  =  2*1  =  2. 

28  gj  =  [2*(0)+l]  [2*(0)+l|  =  1*1  =  1. 

27  gj  =  [2*(0)+l]  l2*(5/2)+!J  =  1*6  =  6. 


ch4 

2 1  1  Here  the  symmetric  hydrogens  give  rise  to  coupling  with 
rotational  wavefunctions  for  the  ratio  of  5:3:2  for  A:F:E 
states.  The  central  carbon  is  what  contributes  to  gj  to 
give  gj  =  [2*(0)+l  |=1. 

3  1  1  Same  as  above  for  coupling  and  gj  =  [2*(l/2)+l]  =  2. 

no2 

446  The  symmetric  oxygens  couple  with  the  rotational 
wavefunctions  to  give  1:0  (even.odd)  ratio.  The  central 
nitrogen  gives  gj  =  [2*(1)+1]  =  3. 
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TEMPERATURE  DEPENDENCE  OF  N2 - BROADENED  HALFWIDTHS  OF  H20 

ROBERT  R.  GAMACHE  and  LAURENCE  S.  ROTHMAN 

The  temperature  dependence  of  N2 -broadened  halfwidths  of 
water  was  studied  for  rotational  vibrational  transitions  of 
atmospheric  importance.  The  halfwidths  were  calculated  using 
cutoff  free  theory  with  Anderson- Tsao-Curnutte  interruption 
functions  and  trajectories  correct  to  first-order-in -time  . 
The  calculations  were  carried  out  for  pure  rotational 
transitions  and  for  v 2  transitions.  The  temperature 
dependence  was  calculated  assuming  the  exponential  form 

7(T„>-7<T)  {|-  }"  • 

The  exponent  n  has  been  calculated  for  different 
rovibrational  transitions  and  the  calculations  have  been 
performed  over  the  temperature  range  200  K  to  1000  K  to  be 
applicable  to  both  atmospheric  studies  as  well  as  hot 
sources.  The  trends  in  n,  such  as  vibrational  dependence,  J 
dependence,  etc.  ,  are  examined  and  recommended  values  are 
given . 
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Temperature  Dependence  Line  Widths  and  Shifts 


Robert  R.  Gamache 

University  of  Lowell,  Center  for  Atmospheric  Research 
450  Aiken  Street,  Lowell,  MA  01854 
and 

Laurence  S.  Rothman 
Air  Force  Geophysics  Laboratory 
Hanscom  AFB ,  MA  01731 


The  temperature  dependence  of  the  pressure -broadened  line  widths  and 
the  pressure  induced  line  shift  are  investigated.  Line  width  and  shift 
calculations  were  done  for  two  major  constituents  of  the  terrestrial 
atmosphere,  ozone  and  water  vapor.  For  ozone  some  126  ro-vibrational 
transitions  belonging  to  either  the  pure  rotation  or  the  i/3  band  were 
studied  via  QFT-ID*  calculations  over  a  temperature  range  of  200-500  K.  H-,0 
has  been  studied*  over  a  temperature  range  of  200-1000  K  to  yield  results 
applicable  to  atmospheric  studies  as  well  as  hot  sources.  Forty-seven 
rotational  transitions  important  to  atmospheric  studies  were  considered  and 
the  calculations  were  done  for  both  the  pure  rotation  and  u2  bands  of  water. 
The  halfwicths  were  calculated  usir.b  a  cutoff  free  theory  with  Anderson- 
Tsao-Curnutte  (h-.TC)  interruption  functions  and  trajectories  correct  to  first 
order  in  time.*  In  both  studies  the  transitions  span  a  wide  range  of  J"  (1 
to  1-  for  H-_ Z  and  1-35  for  03)  with  the  full  manifold  of  Ka"  being  sampled. 


The  calculations  were  done  with  nitrogen  as  the  perturbing  gas  and  the 
calculations  included  the  dipole(K~0  or  03 )  - quadrupole  (N*2 )  ana  the 
cuadrupole  H;0  or  C3 )-  quad-,  upule  (l<2 )  interactions.  The  vibrational 
dependence  was  explicitly  considered  in  all  the  calculations  and  a  Lennard- 
Icr.es  potential  was  used  to  give  the  first-order- in- time  trajectories.  The 
molecular  constants  used  in  the  calculations  ’  *  were  the  best  available  in 
the  literature.  In  these  calculations  all  terms  were  included  (real  and 
imaginary, >  and  no  adjustment  of  the  molecular  parameters  parameters  was 


2y  evaluating  the  halfwidth  at  a  series  of  temperatures,  the  variation 
temperature^  can  be  determined  in  terms  of  a  particular  model.  The 
Tech  in  cm’ ‘/atm  is  expressed  as  the  product  of  density,  velocity,  and 
c  mica 1  cross-section, 


1'  - , 
temperature 


p(~) 


1  oc  i  ty ' v3kT/rru  )  * ,/* )  are  known.  Here  we  assume  the  temperature  dependence 
the  optical  cross-section  follows  a  power  law,  i.e.  cr(T)=-<70  •T"1,  where  aa 
independent  of  temperature.  The  halfwidth  from  Eq .  (1)  is  considered  at 

c  tempera  Cures ,  T.  and  T2 ,  and  the  ratio  is  taken  to  give 


v  (  T  )  •  o  ( T ) 
dependence 


(1) 


o: 


the  density(n0 • 273 . 15/T)  and 


r.e 


l'Ti> 

7(T2) 
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Secring  -n--i/2 
7(Tx) 


m  produces  the  usual  formula 
~  >  -n 


7(T2)  • 


(2) 


The  temperature  dependence  of  the  halfvidch  is  contained  in  the  value  of  one 
exponent  n.  This  model  also  gives  the  temperature  dependence,  a,  of  the 
optical  cross  -  section  which  can  be  useful  for  other  studies. 

A  similar  power  law  form  was  tried  for  the  line  shift, 


A(T) 


,b 
a  a 


(*) 


■owever  we  have  found  that  for  some  transitions 
i(T)-ae^>'L,  fits  the  calculated  line  shifts  better. 


an  exponential 


determine  tne  temperature  exponent,  the  halfwidth 


or 


sr.irt 


is 

a”-luated  at  a  number  of  temperatures  and  the  results  leas t -  squares  fitted 
the  particular  model.  For  example,  for  the  halfwidth  a  plot  (linear 
the  set  of  points  ln(T2/Tl)  vs.  ln(7 (Tx ) /t(T; ) )  yields  a 
straight  line  of  slope-n.  Furthermore,  the  correlation  coefficient  of  the 

o  s  urr.  6  c 


regression) 


tit  or 


<-w 


'.e  poir 


:o  a 


straight  line  relates  to 


ilia: 


the 


as: 


:amperature  dependence  of  the  optical  cross-section,  i.e.,  the  closer 
:orrelation  coefficient  is  to  li.l  the  more  valid  the  assumption. 


full  description  of 


the  results,  with  tables,  is  presented  in 
references  1  and  2,  here  we  only  summarize  the  results. 

In  figure  1  we  show  a  typical  input  set  of  points  and  the  resulting 
straight  line  from  regression  for  the  pure  rotation  transition  2  2  0- — 3  1  3 

cf  water  vapor.  Ail  the  points  fall  on  the  straight  line  as  seen  in  the 
figure.  This  is  confirmed  by  the  correlation  coefficient  of  the  fit  which 
was  C. 99598  for  this  transition.  For  all  of  the  transitions  studied  here, 
the  correlation  coefficient  for  the  halfwidth  was  equal  to  or.e  to  three 
significant  figures  confirming  the  assumed  temperature  dependence  of  the 
optical  cross-section,  i.e.,  ag*!111. 

In  figures  2  and  3  the  temperature  exponents  for  the  halfwidth  are 

'.umbers  for  H-0  and  03  respectively.  Note 
:ure  that  would  help  one  rredict  n  for  a 

temperature  dependence  of  the  shift  results 
ir.g  the  power  law  and  exponential  law  fits.  Mote  in  figure  1  tr.e 
power  law  gives  the  best  results  whereas  it  is  the  reverse  in  figure  5. 
Although  the  shift  calculations  are  still  in  a  preliminary  stage,  they  are 
useful  to  show  trends.  The  results  indicate  large  variations  in  tr.e 
multiplying  and  exponential  coefficients  in  both  models  with  transition.  It 
does  not  appear  that  these  coefficients  can  be  predicted  with  any  certainty 
nor  does  the  use  of  an  average  value  seem  warranted.  There  remains  much  to 
be  done  on  this  subject. 


c  vs.  rotations* 

qua 

appears  tc  be 

no 

uiar  transition. 

igures  4  and  5  s 

how 

ssum.ir.a  the  tower 

law 

mmarv 


Temperature  dependence  of  the  halfwidth 


1  58 


For  both  vater  vapor  and  ozone  it  appears  that  che  temperature 
dependence  of  the  halfwidth  cannot  be  expressed  by  a  simple  relation 
dependent  on  J"  or  any  ocher  quantum  number.  There  is  some  vibrational 
dependence,  but  the  variation  with  rotational  transition  overshadows  this. 
The  fluctuation  in  the  results  is  much  greater  for  H-,0  when  compared  with 
the  results  for  03 . 

The  calculations  give  an  average  temperature  exponent  of  n-0.68  and 
r.- 0.76  for  air -broadening  of  water  vapor  lines  and  ozone  lines  respectively. 

It  is  suggested  that  whenever  possible  temperature  exponents  calculated 
for  specific  ro -vibrational  transitions  be  used. 


Temperature  dependence  of  the  line  shifts 

Although  still  in  a  preliminary  stage,  the  results  of  the  shift 
calculations  indicate  the  temperature  dependence  is  best  described  by  a 
power  law  for  seme  lines  and  by  an  exponential  law  for  other  lines.  There 
seems  to  be  much  greater  variability  of  the  line  shift  with  temperature  as 
compared  to  the  halfwidth.  The  dependence  varies  greatly  from  line  to  line. 
The  results  suggest  that  this  may  be  important  in  certain  planetary 
atmospheres  or  hot  plumes. 
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TEMPERATURE  DEPENDENCE  OF  LINE  WIDTHS  AND  SHIFTS 
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The  temperature  dependence  of  the  pressure  broadened 
line  width  and  the  temperature  dependence  of  the  pressure 
induced  line  shift  are  investigated.  The  model  considers  the 
optical  cross-section  to  vary  with  temperature  as  o(T)  - 
Oq* Tm  where  oq  is  independent  of  temperature.  From  this 
model  the  temperature  dependence  of  the  half-width  is  given 


7(t>  -  7<V  {~}  " 

and  the  temperature  dependence  of  the  line  shift  is 


S ( T )  ,T  . -n’ 


Ai/(T)  -  hu(TQ) 


s7 v  It-J 


7  n  expression  (2),  S  ( T )  is  the  sign  of  the  line  shift  at 
temperature  T  and  is  necessary  in  order  to  avoid  mathematical 
discontinuities.  The  temperature  dependence  of  the  line 


w i d  t h  and  shift  is  contained  in  the  exponents  n  and  n' 
respect  ive ly . 

The  temperature  dependence  of  ^-broadened  halfwidths 
has  been  studied  over  a  temperature  range  of  200-500°K  for 
ozone  and  from  200-  1000°K  for  water  vapor.  For  O3  ,  126  ro- 
vib rational  transitions  were  studied  and  both  the  pure 
rotation  and  1/3  transitions  were  included.  For  water  vapor, 
some  hundred  transitions  were  studied  for  pure  rotation  and 
for  1^2  transitions.  Tne  transitions  were  selected  for  their 
importance  in  atmospheric  applications,  to  consider  a  wide 
range  of  the  rotational  quantum  number  J",  and  the  full 
manifold  of  Ka"  for  particular  J"  values.  The  results  yield 
average  temperature  exponents  of  ozone  and  water  vapor  of  n  - 
0.77  and  n  -  0.68,  respectively.  The  results  for  both 
molecules  are  compared  here. 

The  temperature  dependence  of  the  line  shift  for  N  2  - 
broadening  of  O3  are  also  presented.  The  calculated  line 
shifts  are  somewhat  in  question,  however,  only  the  trend  with 
temperature  was  studied.  The  values  of  n'  presented  show 
large  fluctuations.  The  many  problems  observed  are  discussed 
in  this  work. 
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SIMPLIFIED  EXPRESSION  FOR  THE  TOTAL  INTERNAL  PARTITION  SUM 
AS  A  FUNCTION  OF  TEMPERATURE 


ROBERT  R.  GAMACHE ,  ROBERT  HAWKINS,  and  LAURENCE  S.  ROTHMAN 

The  total  internal  partition  sums  have  been  calculated  for  the 
molecules  and  their  isotopic  variants  on  the  HITRAN  database.  The 
calculations  employed  near  complete  sets  of  energy  levels  for  all 
relevant  rotation-vibration  states  generated  from  the  molecular 
constants.  All  degeneracies  were  accounted  for  in  the  calculations. 
The  calculations  were  performed  over  a  temperature  range  of  75K  to 
400K.  The  accuracies  of  the  resulting  partition  sums  will  be  shown. 
For  generalized  applications,  the  internal  partition  sums  have  been 
fitted  to  the  expression 

Q ( T )  =  a  +  b’T  +  c • T 2  +  d*T3, 

where  a,  b,  c  and  d  are  constants  adjusted  by  means  of  a  least-squares 
minimax  procedure.  The  maximum  deviation  of  the  fit  is  on  the  order 
of  a  few  tenths  of  a  percent. 
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The  total  internal  partition  sums  have  been  calculated  for  the 
molecules  and  their  isotopic  variants  on  the  HITRAN  database.  The 
calculations  employed  near  complete  sets  of  energy  levels  for  all  relevant 
rotation-vibration  states  generated  from  the  molecular  constants.  All 
degeneracies  were  accounted  for  in  the  calculations.  The  calculations  were 
performed  over  a  temperature  range  of  75K  to  400K.  The  accuracies  of  the 
resulting  partition  sums  will  be  shown.  For  generalized  applications,  the 
internal  partition  sums  have  beenfit  to  the  a  simple  four  parameter  Taylor 
series  expansion.  The  maximum  deviation  of  the  fit  is  on  the  order  of  a 
few  tenths  of  a  percent.  This  scheme  will  be  used  on  future  editions  of  the 
HITRAN  database  to  enable  users  to  simply  make  temperature  corrections  to 
the  line  intensities  and  retain  information  about  isotopes  as  well. 
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Temperature  Dependence  of  N2-Broadened  Halfwidths  of  Ozone 

R(  IBI  K  I  R  G\m*uhi 

’  ’  ■  lv*;, Rc'curii  flu  l  ri\cr\ih  c'  J 11  l.'M'x*  .X lrcc: 

\Ltisii  "nuJI*  ;//,S 

I  he  t-.-mp-.-rdUirc  Jt-pciufenee  of  N-hrutiJcned  hdl fwidlbv  ol  n/v'iu  tux  tvvri  studied  over  a 
icrtipcrjiure  rungs'  ol  2( w " 1  K  Xsxunnng  the  lorm  ■>(  7  i  -  •**  /,,«  I  r.  the  exponent  »  has 
lx-er,  oeiet  mined  tr*>m  OPT -ID  eaieulations  tor  |2(>  ro-whraimnal  transitions  ol  cuone  selected 
to  eons Mer  u  mde  range  o!  the  rotational  quantum  number  ./"  II  to  ?5i.  the  lull  manifold  of 
R'l  lor  particular  ./".  salues  and  transitions  important  to  stratospheric  studies  1  he  results  give  an 
average  temperature  exponent  of  n  -  0  In:  nitrogen-broadening  The  temperature  exponents 
are  transition  dependent  hut  the  fluctuation  about  the  average  value  is  small  I The  value 
of  r  as  a  function  of  ./'  and  A  ",  is  investigated  C  ompanson  of  the  present  results  xxuh  expenment 
and  other  theoretical  calculations  shows  good  agreement  Based  on  the  results  an  average  tem¬ 
perature  exponent  lor  air-hroudentng  ol  *'  7  li  :0  is  suggested  c  uss  AcaaemiL  grew.  in. 


I  INTRODUCTION 

The  use  of  balloon-borne  and  ground-based  spectrometers  to  study  and  monitor 
the  atmosphere  has  greatly  increased  our  knowledge  of  mans  regions  of  the  atmosphere. 
It  is  well  recognized  (/.  .1  that  tn  order  to  reduce  the  data  from  these  methods  an 
accurate  knowledge  of  the  broadening  coefficients  for  ro-'ibrational  transitions  is  nec¬ 
essary  For  atmospheric  applications  the  temperature  dependence  of  the  halfwidth 
must  be  understood  as  well  as  the  halfwidth  at  a  given  temperature  for  all  species 
under  consideration.  Here  we  are  concerned  with  the  molecule  ozone  which  has  a 
xerx  neb  spectrum  and  is  known  toplas  an  imoortant  role  in  the  temperature  regulation 
of  the  lower  atmosphere,  air  chemical  cycles,  and  climate  and  which  must  be  well 
determined  in  order  to  extract  information  on  trace  atmospheric  constituents  in  the 
tar  IR  (.'). 

A  fair  amount  of  effort  has  been  invested  in  the  determination  of  accurate  haifwidths 
for  ro-vibralional  transitions  of  o/.one  perturbed  by  atmospheric  gases,  i.e..  N;.  CT. 
air.  Ox.  Several  experimental  measurement  efforts  have  been  performed  14-10)  yielding 
hallwidths  for  some  one  to  two  hundred  transitions,  mostly  for  N;-broadening  with 
some  O;-.  air-,  and  self-broadening  results  as  well.  Theoretical  calculations  using  the 
conventional  ATC'  method  (11)  have  been  performed  (12-15).  The  last  two  of  these 
studies  ( 14.  15).  using  the  most  up  to  date  molecular  constants,  produced  results  that 
were  too  low  compared  with  experiment.  Recently  the  QFT-ID  method  (16)  was  used 
to  study  127  transitions  of  ozone  broadened  by  N;.  0:.  and  air:  comparison  with 
expenment  gave  agreement  to  7', 7.  Using  this  method,  halfwidths  have  been  calculated 
1 1  ~)  for  N;-broadening  of  ozone  for  all  the  unique  rotational  transitions  of  ozone  with 

n  oo::-:85:/85  $3.00 
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A  description  and  summary  of  the  latent  edition  "f  the  A!  (11.  MI  !  HAN  meieiular  h ! »r | »t i« »r.  parameter* 
database  are  presented.  This  ness  database  t  Minium*  the  miurmai :■  >n  forthe-rsen  print' i p r*  1  aimo-phtrif 
absorbers  and  iw**nt\ -one  additional  nn'lei  ular  specie*  [»reviou*ls  nmtamed  on  the  AFGL  atmospheric 
absorption  line  parameter  compilation  and  on  tin-  trace  gas  c» anpilat ion.  In  addition  to  updating  the 
parameters  on  earlier  edit  ion*  ot  t  he  ion,  pda  non.  ness  para  met  er*  has  <•  been  added  to  this  edit  ion  *ui  h  u>  t  he 
>elt  -  broadened  haltsvidth.  the  temperature  dependence  oj  the  air  nroadened  halfssidlh.  and  the  transition 
prohabihiv  Tin-  <lalaha*e  contain*  **  entries  between  0  and  l~.9oo  cm"  .  A  I  ‘  ’ l<  I  i<\\  program  is  nose 
furnished  to  allow  rapid  a<te*>  to  t!ie  moiei  ular  transition*  and  lor  the  creation  ol  customized  output  A 
separate  lib  ol  moiec  ular  <  ross  sections  ot  eleven  heavv  molecular  spec  le*.  applicable  tor  qualitative  simula¬ 
tion  oj  transini**ion  and  c  mission  in  the  ;»t  niosph*  re.  ha*  also  been  pros  tiled 


I.  Introduction 

The  high-resolution  transmission  molecular  absorp¬ 
tion  database  ( known  under  the  acronym  HITRAN  i  is 
a  compilation  of  spectroscopic  parameters  from  which 
a  wide  variety  of  computer  simulation  codes  are  able  to 
calculate  and  predict  the  transmission  and  emission  of 
radiation  in  the  atmosphere.  This  database  is  a  prom¬ 
inent  and  long  running  effort  established  by  the  Air 
Force  at  the  Air  Force  Geophysics  Laboratory  ( AFGL) 
in  the  late  1960s  in  response  to  the  requirement  of  a 
detailed  knowledge  of  infrared  transmission  proper¬ 
ties  of  the  atmosphere.  With  the  advent  of  sensitive 
detectors,  rapid  computers,  and  higher  resolution 
spectrometers,  a  large  database  representing  the  dis- 
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cretc  molecular  transitions  that  affect  radiative  propa¬ 
gation  throughout  the  electromagnetic  spectrum  be¬ 
came  a  necessity.  A  wide  range  of  applications  for 
HITRAN  has  evolved  including  detection  of  trace  and 
weakly  absorbing  features  in  the  atmosphere,  atmo¬ 
spheric  modeling  efforts,  laser  transmission  studies, 
remote  sensing,  lidar,  and  a  reference  base  for  funda¬ 
mental  laboratory  spectroscopic  research.  The  HI¬ 
TRAN  database  has  been  periodically  updated  and 
enhanced  since  it  first  became  generally  available.’-4 
The  most  recent  edition  of  the  HITRAN  database  was 
made  available  in  late  1986.  This  latest  version  now 
unites  the  data  on  twenty-eight  molecular  species  with 
bands  covering  regions  from  the  millimeter  through 
visible  portion  of  the  spectrum.  Originally  the  data¬ 
base  contained  for  each  molecular  transition  the  fol¬ 
lowing  basic  parameters:  (1)  resonant  frequency;  (2) 
line  intensity;  (3)  air-broadened  halfwidth;  and  (4) 
lower  state  energy  (as  well  as  unique  quantum  identifi¬ 
cations).  Additional  parameters  have  recently  been 
provided  which  permit  new  capabilities  for  remote 
sensing  in  the  atmosphere  and  capabilities  to  deal  with 
nonlocal  thermodynamic  equilibrium  effects  in  the  up¬ 
per  atmosphere.  The  overall  structure  of  the  database 
has  been  expanded  to  include  files  of  cross-sectional 
data  on  heavy  molecular  species  such  as  the  chlorofluo- 
rocarbons  (CFCs)  and  oxides  of  nitrogen  which  are  not 
yet  amenable  to  line-by-line  representation.  This  has 
added  to  HITRAN  the  capability  of  qualitative  detec¬ 
tion  of  anthropogenic  gases  in  the  window  regions  of 
the  infrared.  Ongoing  research  efforts  will  gradually 
move  some  of  these  data  to  the  main  body  of  the 
database.  The  new  file  struct  ure  of  HITRAN  is  shown 
in  Fig.  1.  New  parameters  have  been  added  to  the 
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Temperature  Dependence  of  NrBroadened  Halfwidths  of  Water 
Vapor:  The  Pure  Rotation  and  v2  Bands 

Robi  r  i  R.  GamaChi: 

inner (  *T»,V-r  '»  |  ,•#»,*.  I'lf  .l;v';  .\/r,v:  /.«*nt,.V  W.;' 'JT/iiniT.G  l/lt.'-f 

AM) 


I.M'Ri  so  S.  Rothman 

\tH  iiUH-H' i‘,‘ ! 

1  Ik-  temperature  dependence  nt  N:-bmadencd  halfwidthx  ol  H;C)  hav  been  studied  over  a 
temperatuie  range  uf  « <-  HKio  K  u>  ucld  results  applicable  to  atmospheric  studies  as  well  as  hoi 
sources  1  ortv-seven  rotational  transitions  important  to  atmospheric  studies  were  considered  and 
the  calculations  were  done  lor  both  the  pure  rotation  and  the  i;  bands  ot  water  The  transitions 
span  a  range  ot  J‘  from  I  to  1-1  with  the  manifold  of  AT  being  sampled  The  hall’vsidths  were 
calculated  using  a  cutoff-free  theors  vsith  Anderson-Tsao-Curnutle  interruption  functions  and 
trajectories  correct  to  firs:  order  in  time  Assuming  the  form  •><  1 1  -■  si  7C. K  I\,:  I )'.  the  exponent 
n  has  been  determined  lor  each  transition  I  he  results  show  strong  fluctuations  in  r,  with  the 
rotational  quantum  numbers  and  a  lesser  dependence  on  the  particular  vibrational  hand  The 
results  give  average  temperature  exponents  for  V  broadening  ol’«(pure  rotation)  =  U.6K4  and 
hii;)  -  li.h'h  Based  on  the  results,  an  average  value  of  tt  =  0.68  is  recommended  and  it  is 
suggested  that  values  calculated  for  specific  rotational  transitions  be  used  when  available.  - 
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I  IN TRODl'C  TION 

The  use  of  spectroscopy,  via  balloon-borne  or  ground-based  spectrometers,  has 
greatly  increased  our  knowledge  of  the  constituents  and  chemistry  of  the  atmosphere. 
The  successful  retrieval  of  stratospheric  information  from  remotely  sensed  data  requires 
that  accurate  spectral  parameters  (e.g..  line  intensities,  line  widths,  and  line  frequencies! 
be  known  t/i.  In  addition  to  these  parameters,  the  temperature  variation  of  the  pa¬ 
rameters  must  be  known,  from  200  to  .TOO  K.  t. -r  atmospheric  studies  and  in  the  range 
of  1000  K  for  hot  sources  This  paper  focuses  on  v-.er  vapor.  Although  a  minor 
constituent  of  the  terrestrial  atmosphere,  water  vapor  has  a  very  rich  spectrum  and  is 
ol  major  importance  for  considering  the  absorption  of  infrared  radiation  in  the  at¬ 
mosphere  1 2).  It  plays  a  prominent  role  in  determining  atmospheric  transmission  to 
solar  or  laser  radiation  and  has  a  profound  effect  on  the  heat  balance  of  the  lower 
atmosphere  i.h  This  species  is  of  primary  concern  in  ihe  characterization  of  hot 
sources  furthermore,  in  order  to  extract  information  about  trace  gas  constituents  of 
the  atmosphere,  the  principle  absorbers  must  be  well  characterized  (•/). 

for  water  vapor.  a  fair  amount  of  effort  has  gone  into  determining  halfwidths  for 
a  number  ol  rowbraiinnal  transitions  both  experimentally  (.'-/G  and  theoretically 
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